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R.  Eugene  Collins,  Pay  E.  Calvert,  and  Humphrey  H.  Hardy 
University  of  Houston,  Houston,  Texas 


Progress  Report,  December  1979 


The  problem  of  developing  a  computer  simulator  for  the  human 
cardiopulmonary  system  has  been  approached  from  two  directions  in 
the  first  phase  of  our  project,  development  of  a  pulmonary  model 
and  development  of  a  circulatory  model,  each  distinct  from  the 
other.  The  second  phase  of  this  effort  addresses  interfacing  of 
these  two  models  into  a  single  functional  model  and  i ncroporation 
of  simulation  of  physiologic  feed-back  and  control  mechanisms  into 
this  model.  This  report  consists  of  three  sections  ana  describes 
the  major  progress  and  results  achieved  in  the  first  phase  of  this 
project. 

Summary  of  Section  A 

Section  A  is  the  text  of  a  paper  being  submitted  for  publicati 
which  concisely  describes  the  development  of  the  equations  for  the 
circulatory  model  as  well  as  the  overall  progress  made  in  using  this 
model  to  examine  the  effects  of  Whole  Body  Acceleration  (WBA)  on  the 
circulatory  system.  Greater  detail  on  the  development  e  :  tn’s  model 
was  given  in  our  last  annual  report. 

Summary  of  Section  B 

Section  B  consists  of  the  text  of  a  paper  being  submitted  for 
publication  which  describes,  in  detail,  our  handling  of  the  pressure 
volume  relationship  desci ibing  the  elastance  of  the  circulatory 
elements.  The  relationship  employed  is  significant  in  that,  though 


2 


it  is  achieved  by  curve  fitting  techniques,  it  not  only  provides  an 
excellent  fit  to  existing  data  on  circulatory  elements,  but  it  is  a 
relationship  that  has  been  employed  in  describing  elastance  in  many 
other  organic  tissues  and  may  represent  a  relationship  of  much 
broader  applicability. 

Summary  of  Section  C 

Section  C  consists  of  the  text  of  R.  E.  Calvert's  doctoral 
thesis  which  includes  a  description  of  the  overall  development  of 
the  pulmonary  model  as  well  as  the  major  progress  made  during  the 
last  year.  An  indepth  examination  of  the  combined  effects  of  WBA 
and  seatback  angle  on  pulmonary  mechanics  is  presented.  A  principal 
finding  of  this  study  being  the  prediction  that  a  seatback  angle  of 
between  60°  and  75°  should  be  optimal  in  relieving  the  adverse  effects 
of  WBA  on  respiration. 

A  more  complete  summary  of  each  section  may  be  found  in  the 
abstracts  preceeoing  each  section. 

Phase  one  of  this  project  is  nearly  complete  and  the  separate 
models  are  functioning  and  pi'oviding  interesting  results  describing 
their  respective  systems.  Phase  two  has  been  initiated  during  the 
past  year.  This  phase  consists  of  interfacing  the  two  models  and 
developing  the  overlaping  feedback  systems  which  will  describe  the 
functioning  of  the  entire  respiratory  system.  Interfacing  of  the 
pulmonary  model  with  the  pulmonary  circulation  portion  of  the  circula¬ 
tory  model  has  been  achieved  and  some  results  obtained.  (See  Section  A) 


SECTION  A 


The  Circulatory  Model 


A  DIGITAL  COMPUTER  MODEL  OF  ThlE  HUMAN 


CIRCULATORY  SYSTEM* 

H.  H.  Hardy,  R.  E.  Collins,  and  R.  E.  Calvert 
The  University  of  Houston,  Houston,  TX  77Q04 


ABSTRACT 

A  digital,  many  compartment  model  of  the  human  circulatory 
system  has  been  developed  which  simulates  pulsatile  olood  flow  and 
gas  transport  and  exchange.  The  model  is  designed  specifically  to 
study  short  term  whole-body  acceleration  (g^)  encountered  in  modern 
aerial  combat  maneuvers  and  incorporates  a  realistic  representation 
of  the  nonlinear  elastic  characteristics  of  circulatory  eleirients  and 
the  related  pressure  dependent  flow  resistance  characteristic  of  these 
elements. 

The  model  has  been  shown  to  properly  simulate  numan  data  for 
passive  breathing  in  a  supine  subject.  The  computed  carbon  dioxide 
and  oxygen  partial  pressures  vary  realistically  around  measured 
average  partial  pressures  for  human  subjects.  The  computed  time 
variation  of  pressures,  circulatory  chamber  volumes  and  flow  rates 
match  corresponding  human  data.  Under  sinusoidal  g  variations,  tne 
model  predicts  realistic  variations  of  volumes,  flows  and  pressures. 


^Supported  in  part  by  the  Air  Force  Office  of  Scientific  Research. 
Grant  Number  AFOSR-76-2905. 


INTRODUCTION  AND  REVIEW  OF  EXISTING  MODELS. 


A  number  of  models  of  the  circulatory  system  in  man  have 

been  developed.  In  the  late  fifties  and  sixties,  great  strides  were 

made  by  many  researchers  in  an  effort  to  model  the  circulatory 

system.  The  analog  computer  was  well  developed  and  provided 

modeling  of  heart  dynamics,  which  could  then  be  used  as  a  diagnostic 

aid.  This  work  reached  a  peak  in  1962  when  Circulatory  Analog 

Computers^^^  was  published  from  the  proceedings  of  the  symposium  on 

the  applications  of  analog  coinputers  to  the  study  of  the  mammalian 

circulatory  system  held  in  the  Netherlands  on  April  19  and  20,  1962. 

Following  this,  study  turned  from  the  circulation  itself  to  the 

regulation  of  the  circulation.  In  1966,  a  conference  was  held  on  the 

regulation  of  the  circulation  and  metabolite  transport  resulting  in  the 

book.  Physical  Bases  of  Circulatory  Transport;  Regulation  and 

(2) 

Exchange.  Since  that  time  research  in  all  three  areas  (arterial 
pulse  propagation,  regulation,  and  gas  exchange)  has  continued 
independently. 

A  summary  of  some  of  the  circulatory  models  that  have  been 
developed  is  given  in  Table  1  (this  model  is  included  in  the  table  for 
comparison).  The  table  indicates  the  author,  date  of  publication,  in 
what  detail  the  systemic  and  pulmonary  systems  have  been  divided,  to 
what  extent  physiologic  data  has  been  linearized,  and  the  studies 
performed  with  the  models.  From  the  table,  it  can  be  seen  that 
existing  models  of  the  circulatory  system  can  be  divided  into  three 
distinct  categories: 


1 


2 


1.  Those  models  containing  a  pumping  heart--these  are  all  analog 

in  nature  and  include  Dick  and  Rideout/^^^  Noordergraaf/^ 

^  ,  (10-12)  ^  (14)  ,  „  (15) 

beneken,  Hwang,  and  Croston. 

2.  Those  models  addressing  the  effect  of  acceleratory  stress— only 

(18) 

two  are  known  to  the  authors:  that  of  Green  and  that  of 
Boyers. 

3.  Those  models  containing  gas  transport  and  diffusion— many  lung 

models  liave  some  form  of  circulation,  usually  this  consists  only 

of  an  influx  of  veneous  blood  and  an  efflux  of  arterial  blood. 

These  are  not  considered.  Those  models  with  a  more  complete 

(21,22)  (19) 

circulatory  system  include  those  of  Saidel,  ’  Grodins. 

w  P  ^20) 
and  cmerv 


No  model,  known  to  the  authors,  fails  wiilNn  any  two  of  these 

categories.  The  model  developed  here  addresses  all  three  categories  and 

therefore  is  a  more  comprehensive  model. 

In  addition,  only  one  of  the  models  represents  flow  resistances  as 

dependent  on  vessel  geometry  as  required  by  the  well-established 

(24) 

Poiseuille-Hagen  relationship.  All  but  two  of  the  models  assume  linear 
pressure-volume  or  pressure-diameter  relationships,  yet  all  in-vivo  data  and 

inodern  in-vitro  data  (holding  vessel  length  constant)  indicate  a  semi- 

,  u  •  .  u-  (27-31) 

logarithmic  relationship 

The  one  model  with  variable  resistance  is  only  an  arterial  model  and 
IS  thereiore  not  useful  to  study  the  entire  circulatory  system.  The 
assumptions  made  in  the  remaining  models  restrict  their  reliability  to  small 
pressure  changes  and  make  them  of  questionable  value  in  studying  g-stress. 
The  model  presented  here  will  incorporate  better  pressure-volume 
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TABLE  1.  Summary  of  Charac ter i s L i os  of  I'ul)  1  i shod  I'loiJels 


and  resistance-volume  relationships  than  any  compleie  circulatory  model 
now  available  in  order  to  make  possible  a  more  accurate  study  of  g-stress. 
In  addition,  it  will  include  the  exchange  and  transport  of  gases  with  a 
pumping  heart. 

The  model  presented  here  also  includes  the  following  anatomical 
details  to  simulate  further  the  human  circulatory  system: 

1.  Venous  return  volume  has  been  separated  from  the  capillary 
volume  so  that  venous  storage  without  significant  and  CO^ 
exchanges  with  tissue  is  possible. 

2.  Each  lung  and  body  compartment  has  venous  return  to  provide 
for  the  gravitational  effect  of  different  body  orientations  and 
venous  storage. 

3.  The  leg  return  compartment  and  valves  simulate  control  of 
circulation,  and  muscle  action  on  circulation,  within  the  legs. 

4.  The  chest  cavity  compartment  allows  changes  in  blood  pressure 
and  flow  in  response  to  changes  of  pleural  and  alveolar  pressure 
due  to  breathing. 

5.  .Multiple  vascular  compartments  in  the  lung  will  simulate 
distribution  of  perfusion  and  gas  exchange  with  air  in  the  lung. 

6.  External  pressure,  P',  and  resistances  at  each  compartment  are 
variable  to  simulate  short  term  neural  control,  g  suit  effects, 
and  breathing  maneuvers. 

7.  Each  portion  of  the  body  has  a  separate  control  of  resistance 
and  compliance  allowing  simulation  of  vasoconstriction  and 
dilation  in  the  skeletal  muscles  while  allowing  no  vasoconstric¬ 
tion  ir.  cranial  arteries  and  neither  vasoconstriction  nor  vasodil- 
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ation  in  the  lungs. 

8.  A  detailed  simulation  of  heart  function  which  permits  studies  of 
such  abnormalities  as  a  defective  heart  valve,  high  blood 
pressure,  or  enlarged  heart. 

9.  The  variable  resistance  of  the  distensible  tubes  of  the  vascular 
system  is  adequately  modeled. 

10.  The  dependence  of  compartment  pressures  on  gravity  and  whole 
body  acceleration  is  included  for  ail  body  orientations  (e.g., 
horizontal,  reclining,  legs  up,  head  down,  etc.) 


BLOOD  FLOW  EQUATIONS 

Here  we  address  the  mechanics  of  blood  circulation  in  the  model. 
The  discussion  of  the  blood  flow  equations  is  divided  into  three  parts.  The 
first  part  treats  the  axial  motion  of  the  blood  through  the  blood  vessels,  the 
second  part  treats  the  transverse  motion  of  blood,  vessel  wall,  and 
surrounding  tissue,  and  the  third  part  treats  the  vascular  pressure-volume 
and  resistance  relationships  used  in  the  model. 

Throughout  the  model,  the  basic  assumption  is  that  many  branching 
blood  vessels  may  be  approximated  by  a  single  vessel  of  uniform  size.  Due 
to  the  enormous  geometric  complexity  and  number  of  blood  vessels  in  the 
body,  this  assumption  is  essential  if  a  model  is  to  be  developed.  The  single 
vessel  in  turn  is  mathematically  equivalent  to  a  capacitive  element 
(henceforth  referred  to  as  a  vessel  or  model  chamber)  and  a  resistive 
element  (henceforth  referred  to  as  a  model  segment).  The  combination  of 


the  two  will  be  referred  to  as  a  model  element. 


These  model  elements  are  linked  together  to  form  the  entire 
circulatory  model  (Figures  1  and  2).  Only  three  types  of  capacitive 
elements  are  employed  in  the  entire  model.  Type  a)  consists  of  a  blood 
volume  element  able  to  exchange  gases  with  an  alveolar  space  in  the  lungs 
(e.g.:  Chamber  5,  Fig.  2),  b)  is  a  blood  vessel  able  to  exchange  gas  with  a 
tissue  volume  element  (e.g.:  Chamber  15,  Fig.  !),  and  c)  is  a  blood  volume 
element  having  no  exchange  of  gas  with  exterior  elements  (e.g.:  Chamber 
1,  Fig.  1).  Types  a)  and  b)  will  be  referred  to  as  the  capillary  elements. 

Valves  are  indicated  in  Figure  1  within  the  heart  and  leg  chambers. 
The  flow  of  blood  through  the  leg,  mitral  and  tricuspid  valves  is  modeled  by 
the  specifications  such  as: 


(1) 


^20  ^19’  '^19  ^  ^  valve  open 


if  P^Q  valve  closed 


where  v^^  is  the  volumetric  flow  rate  of  blood  in  ml/sec  from  chamber  19 

into  20,  with  pressures  P,^  and  P.^  respectivelv  m  mm  11  g  (sec  Figure  1. 

19  20 

legs  and  leg  return).  This  allows  only  one-way  flow  through  the  leg  valves. 

The  aortic  and  pulmonary  valves  also  allow  only  one- wav  flow,  but 
are  modeled  in  a  slightly  different  way  due  to  Inclusion  of  the  inertia  of  the 
blood.  If  the  aortic  or  pulmonary  valve  is  closed,  the  same  conditions  apply 
as  in  Equation  1.  If  the  valve  is  open,  then  the  volumetric  flow  rate  is 
calculated  from  the  equations  3,  5,  and  6,  and  its  sign  is  checked;  if 
positive,  calculation  continues;  if  negative,  the  volumetric  flow-  rate  is  set 
to  zero,  i.e.,  the  valve  is  closed. 


>- 


o  oc. 
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Figure  2.  I'ulmonury  ciixjulatioii  ot  Lht.*  iiiodt: 
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AXIAL  MOTION 


The  equations  for  axial  tnotior.  of  the  blood  through  a  circulatory 

element  are  derived  from  momentum  conservation,  written  in  terms  of 

forces  as  source  terms.  The  forces  acting  consist  of  the  pressure  force  of 

the  fluid  at  the  inlet,  and  the  outlet,  -^^2^2'  fot'ce  of  gravity, 

pg(h  -h  )A,  and  the  viscous  drag  force  at  the  vessel  walls,  -f  A  with  A 
12  s  s 

the  surface  area  of  the  element.  Here  A  is  the  total  surface  area  of  the 

s 

vessel;  A  is  the  average  cross  sectional  area  of  the  vessel;  is  the 

distance  in  cm  separating  the  inlet  and  outlet  ends  of  the  vessel  along  the 
direction  of  the  gravitational  force;  and  f  is  the  total  drag  force  per  unit 
area  of  vessel  wall.  The  algebraic  sum  of  these  forces  is.  according  to 
Newton's  second  law,  a  source  of  momentum  in  a  conservation  equation  for 
the  axial  momentum  of  the  blood  in  the  vessel.  In  terms  of  the  mean 


volumetric  flow  rate  through  the  vessel,  q,  inlet  and  outlet  rates,  q^,  q^  or 
velocities,  v^,  v^,  and  the  length  of  the  vessel,  L,  the  equation  of  motion  of 
the  blood  is  then 

(2)  *  =  -’'i'll  -  '^''2''2  '  '’Al  -  ‘’2''2  *  ■  "2’  ■  ‘''s 

Within  each  model  element  the  axial  flow  segment  is  assumed  to  have 

a  uniform  cross  sectional  area  (i.e.,  A^  -  A^  =  A)  and  constant  vessel 

length,  L.  The  contribution  of  drag  forces,  fA  ,  can  be  approximated  by 
2 

(R'q  +  Rq)  A  where  R'  and  R  are  functions  of  the  radius  of  the 
(25) 

vessel  .  Additional  simplifying  assumptions  are  made  for  each  of  the 
model  elements.  These  assumptions  are  best  exhibited  if  the  above 
equation  is  written  as; 


(pL)dq 
A  dt 


2 

P  Vj  +  pv^ 


-  Pq  =  ?j  -  ^2 


Pg(h,  - 


(3) 


2 
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This  is  the  equation  of  axial  motion  of  blood  in  the  model  segments 
connecting  left  ventricle  and  aorta  and  connecting  right  ventricle  and 
pulmonary  artery  where  the  cross  sectional  area  is  assumed  constant. 

In  all  other  model  segments,  the  inertial  terms  in  Equation  (3)  are 


considered  negligible  compared  to  the  drag  force  term,  i.e., 
(4)  1  2 


2  2  2 

(pL)dq  -  Pv  +  Pv  «  R’q  +  Rq 
A  dt 


This  assumption  is  justifiable  since  the  rate  of  change  of  the  volumetric 

fluid  velocity  is  small  compared  to  the  rate  of  volumetric  flow  under 

normal  conditions.  With  rapidly  changing  g,  this  assumption  must  be 

reassessed,  but  enough  information  is  not  yet  available  for  a  conclusion  at 

this  time.  It  is  also  assumed  that  in  these  segments,  the  quadratic  term 
2 

R'q  is  negligible  compared  to  the  viscous  term  Rq.  This  is  reasonable 

because  flow  rates  are  appreciable  only  between  the  ventricles  and  arteries 

(26) 

at  the  peak  of  systole.  These  assumptions  yield  the  equation  for  axial 
flow  in  these  segments  as, 

(5) 

In  addition  to  the  above  equations,  conservation  of  mass  for  an 


Pi  -  P2  ^  Pg(hj  - ^  Rq 


incompressible  fluid  yields  an  equation  applicable  to  all  of  the  model 

dV 

chambers  for  the  change  of  chamber  volume  ~  in  terms  of  the  flow  into 

dt 

the  chamber,  v^,  and  the  flow  out  of  the  chamber  v^: 

1  dV 


(6) 


A  dt 


transverse  motion 


The  equations  for  transverse  motion  of  the  vessel  wall  and  surround¬ 
ing  tissue  are  also  derived  from  Newton's  second  law.  The  assumption  is 
made  that  there  is  no  pressure  gradient  tranverse  to  the  vessel  axis  within 
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each  of  the  vessels.  An  additional  assumption  of  instantaneous  equilibrium 

for  the  transverse  motion  of  the  vessel  walls  is  made.  Thus  the  blood 

pressure  inside  the  vessel,  P,  is  equal  to  the  external  pressure  on  the  vessel 

and  surrounding  tissue,  P',  plus  the  elastance  pressure  from  the  muscular 

tissue  surrounding  the  vessel,  P,  plus  a  resistive  term  consisting  of  a 

constant,  R,  times  the  rate  of  change  of  the  volume  of  the  vessel.  Tiius 

^  -  -  dV 

(7)  P.P..P.R- 

This  is  the  equation  used  in  the  model  for  the  transverse  motion  of  the 
vessel  wall  and  surrounding  tissue  in  the  arterial  chambers.  The  assumption 
of  instantaneous  equilibrium  was  made  after  an  in  depth  computer  study 
which  indicated  that  this  assumption  affected  only  the  dicrotic  notch  and 
had  negligible  effects  (^'5  mm  Hg)  on  the  rest  of  the  arterial  pulse  in 
modeling  a  passive,  supine  position. 

In  order  to  simplify  calculations  in  the  other  chambers  of  the  model, 
the  viscous  term  was  considered  negligible.  The  justification  for  the 
removal  of  this  term,  as  with  the  assumption  of  instantaneous  equilibrium, 
can  only  be  established  by  operating  the  model  under  various  conditions  and 
comparing  the  results  to  physiologic  data.  Other  models  in  which  the 
viscous  term  was  included  have  demonstrated  that  R  is  very  small,  but 
final  acceptance  of  this  approximation  will  await  many  computer  experi¬ 
ments  on  the  model.  With  these  assumptions,  the  transverse  motion  of  all 
other  chambers  is  represented  as 

(8) 

P  =  P'  +  P  . 
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THE  PRESSURE-VOLUME  RELATIONSHIP  FOR  CIRCULATORY  ELEMENTS 


As  was  mentioned  in  the  section  on  axial  motion,  the  flow  resistance 
factors,  R  and  R',  are  dependent  on  the  vessel  radius,  r.  In  the  chambers 
including  the  non-linear  term  in  q,  however,  the  cross  sectional  area  of  the 
vessels  and  therefore  the  radius  of  the  vessels  is  assumed  constant.  Thus 
R'  is  a  constant  for  the  chambers  in  the  model  in  which  it  is  not  zero  . 


The  dependence  of  R  on  the  vessel  radius,  however,  must  be  included 
throughout  the  model.  This  is  done  using  the  Poiseuille-  Hagen  formula 


r 


with  A'  a  constant.  Assuming  cylindrical  vessel  geometry,  this  formula 


can  be  written  in  terms  of  the  vessel  volume, 

R 

(10)  -  ^2 


with  A  a  new  constant,  since  vessel  length  is  assumed  fixed. 

Since  the  volume  of  the  vessel  is  a  nonlinear  function  of  the 
transmurai  pressure,  an  accurate  P-V  relationship  is  essential  in  setting 
flow  resistances  and  including  g-effects  in  the  model.  Using  in  vivo  data 
from  the  literature  for  dogs  and  man,  it  was  determined  that  the  P-V  data 
for  circulatory  elements  could  be  fitted  more  accurately  by  an  equation  of 


the  form 


-WP  -P  ) 
V-VQ=T(l-e  ^  ^0) 


than  by  the  relationship  of  Bergel  or  variations  of  it  presented  by  other 
(27-31) 

investigators.  Here  (P^  ,V^)  corresponds  to  any  arbitrarily 

selected  Transmural  pressure- volume  point  on  the  curve,  k  is  a  constant 


characterizing  the  elastic  property  of  the  tissue  associated  with  the 
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equilibrium  is  also  assumed  in  the  tissue  chambers  and  a  uniform  gradient  is 
assumed  within  the  vessel  walls  through  which  the  gases  diffuse.  This  is 
not  the  case  in  the  real  system  of  course  and  some  studies  have  addressed 

K, 

this  problem. 

Using  these  assumptions,  and  the  conservation  of  mass  of  each  gas, 

the  rate  of  change  of  the  number  of  moles  of  each  gas  within  a  vessel 

chamber,  dN  /dt,  is  equal  to  the  number  of  moles  transferred  into  the 
8 

chamber  per  second  by  all  the  blood  entering  the  chamber  (Cq)  ,  minus 
the  number  of  moles  transferred  out  of  the  chamber  per  second  by  all  the 
blood  leaving  the  chamber,  {Cq)°^^,  minus  the  rate  of  diffusion  from  the 
chamber  into  the  tissue,  D(P-P): 


(12) 


dN 

^=[Cq]‘'’-[Cq]““'-D(P-P) 


where  P  is  the  partial  pressure  of  the  gas  in  the  tissue  or  alveoli  (in 
mm  Hg)  adjacent  to  the  capillary  chamber,  and  F  is  the  partial  pressure  of 
the  gas  ip.  the  vessel  chamber.  D  can  be  related  to  the  diffusion  constant, 
0 ,  as  follows; 

D  - 

(13)  ^  '  I 


where  Y  is  the  ratio  of  concentration  to  partial  pressure  within  the  tissue, 
A'  is  the  cross  sectional  area  through  which  diffusion  occurs,  and  I  is  the 
distance  over  which  the  diffusion  occurs.  The  volumetric  blood  flow  into 
and  out  of  every  chamber,  q,  is  provided  by  the  portion  of  the  model 
discussed  in  the  previous  section. 

The  model  has  been  designed  to  model  the  transport  of  any  number  or 
kinds  of  gases  transported  by  the  blood.  Presently  however,  only  and 
CO^  are  included.  Although  there  are  a  number  of  variables  that  affect 


circulatory  element  and  t  is  a  constant  proportional  to  the  tissue  volume 
associated  with  the  circulatory  element. 

Throughout  the  model  the  resistance  factor,  R,  is  a  nonlinear  function 
of  the  vessel  volume  (Equation  10).  This  is  in  turn  a  function  of  the  elastic 
properties  of  the  vessel  wall  and  the  transmural  pressure  across  the  vessel 
wall  as  modeled  by  Equation  11. 

The  blood  flow  equations  are  solved  numerically  with  the  gas 
transport  equations  as  described  in  the  next  section. 

GAS  TRANSPORT  EQUATIONS 

As  was  illustrated  in  Figures  2  and  3  the  following  assumptions  are 
mtroduced  in  order  to  simplify  the  gas  transport  portion  of  the  model: 
Gases  enter  or  leave  the  circulatory  system  only  by  diffusion  from  the 
pulmonary  capillaries  to  the  alveoli,  which  lead  to  the  outside  air  through 
the  airways  of  the  lung,  or  by  diffusion  from  the  capillary  beds  into  body 
tissues  where  a  loss  of  uxygen  and  a  gain  of  carbon  dioxide  occurs  due  to 
tissue  metabolism.  Throughout  the  rest  of  the  circulatory  system  the  gases 
arc  merely  transported,  the  total  number  of  moles  of  each  gas  being 
conserved.  Within  each  chamber  gas  may  be  gained  or  lost  to  other 
chambers  by  two  processes;  blood  flow  or  diffusion.  Within  all  but  the 
capillary  chambers,  the  mass  of  each  gas  may  be  increased  or  decreased 
only  by  the  influx  or  efflux  of  blood  containing  that  gas. 

Perfect  mixing  in  every  blood  and  tissue  chamber  of  this  model  is  also 
assumed.  As  a  consequence,  the  concentration  of  each  gas  in  whole  blood 
leaving  a  chamber  is  always  equal  to  the  concentration  of  that  gas  in  the 
chamber.  Instantaneous  chemical  equilibrium  is  assumed  between  gases  in 
solution  in  plasma  and  bound  to  hemoglobin  within  the  blood.  Instantaneous 
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the  number  of  moles  of  and  CO^  in  the  blood,  so  far  in  the  model  the 

number  of  moles  of  each  of  these  gases  is  a  function  of  only  the  partial 

pressure  of  both  gases  and  the  pH;  the  pH  is  in  turn  considered  a  function 

of  the  partial  pressure  of  both  gases. 

The  number  of  moles  of  in  each  vessel  chamber  is  calculated  as 

p 

the  sum  of  the  number  of  moles  of  dissolved  in  the  plasma,  N  ,  plus 

h  2 

the  number  of  moles  of  O  bound  to  the  hemoglobin,  N  .  The  number  of 
moles  of  dissolved  in  the  plasma  is  directly  proportional  to  the  partial 
pressure  of  oxygen,  ,  and  the  volume  of  plasma,  (l-H)V,  where  H  is  the 
hematacrit  and  V  is  the  volume  of  whole  blood: 

(14) 

Here  a_  is  the  solubility  factor. 

°2 

The  number  of  moles  of  bound  to  the  hemoglobin  in  each 

chamber  ,N^  ,  is  expressed  in  terms  of  the  fractional  saturation  of 

hemoglobin,  5^  .  Since  totally  saturated  blood  contains  .201  volume 
^2 

percent  of  O  (8.98  moles  of  O  per  ml  I'emcglcocr. ,  '.-.’c  Piave 

(.4)  -O,  =  ^ 

2  2 

where  N  =  number  of  moles  of  O  in  a  chamber  of  blood  volume  V  and 

g  =  8.98  moles  O^/ml  heiroglcbiri 

S  =  fraction  of  total  0_  saturation  of  hemoglobin. 

2  ^ 

This  gives  for  the  total  number  of  moles  of  in  each  chamber,  N  , 

(16)  ^ 

2  2  2  2 


and  for  the  concentration  of  O  , 

2 


6 


(17) 


°2  °2  ^2  ^2 


Equation  16  becomes,  for  O^, 


(18) 


dP  dP^ 

°2  =°2 


o. 


where  a  =  a  (l-H)V  +  nVH — r~ 

°  O, 


as 


o. 


ap 


HV 


CO. 


and 


(19) 


in  out  P  -  ,  dV 

^  -  [C,  q]  -D.  [P^  -  P^  ]  -  — 

(34) 

The  expression  to  be  used  for  S  ^  is  the  one  given  bv  Gomez  .  In  the 

O  °  ' 

2 

process  of  this  calculation,  an  equation  relating  pH  to  the  partial  pressure 
01  and  CO^  is  required  and  the  one  used  by  Kelman^^^^  and  '\)Cest^^^' 
was  chosen. 

In  the  case  of  carbon  dioxide  the  gas  content  of  blood  was  not 
separated  into  plasma  content  and  hemoglobin  content  as  was  done  for 


oxygen,  instead  the  number  of  moles,  >  in  a  chamber  is  written  in 

terms  of  the  partial  pressure  of  carbon  dioxide,  P  .  the  volume  of  whole 

CO^ 

blood.  V,  and  a  saturation  function  for  whole  blood  that  was  adapted  from 
Kelman^^^^  and  West^^^\  : 


(20) 


N  ^  S  P  V 

co^  co^  co^ 


The  molecular  carbon  dioxide  concentration  is  then 


C  =  S  P 

CO^  co^  co^ 


(21) 


17 


and  Equation  20  becomes 


where  c  = 


=  [C_  q]‘"-  [C  o]°'"  -C_ 


CO^  dt 


Equations  17  and  23  are  solved  simultaneously  yielding 

dP  dK  -  bK 

Ji2.  ^  f  °J 

(24)  dt  ad  -  be 

dP^^  aK^^  -cK^ 

_ff2 _ ^ 

(25)  dt  ad  -  be 

Equations  24  and  25  are  used  to  model  the  gas  transport  into  and  out 
of  a  model  chamber  and  are  coupled  to  the  blood  transport  equations 
discussed  in  the  previous  section  and  tissue  chamber  equations  discussed 
below. 

The  amount  of  oxygen  and  carbon  dioxide  in  each  of  the  tissue 
chambers  is  altered  by  the  exchange  of  each  gas  with  capillary  chambers  by 
diffusion  and  by  the  disappearance  of  oxygen  and  the  appearance  of  carbon 
dioxide  due  to  tissue  metabolism.  The  tissue  chamber  equ  tions  are 


therefore 


•  +  D^  (P  -  P  )  -  VI 
°2  °2  °2 


(27) 


dN 


CO, 


dt 


^co  ^  ^ 

2  2  2  2 


where  and  are  the  number  of  moles  of  oxygen  and  carbon 

dioxide;  and  and  rates  of  disappearance  of  and 

^2  ^2 

appearance  of  CO^  ,  respectively. 

In  the  tissues,  the  number  of  moles  of  each  gas  is  assumed  to  be 
proportional  to  the  partial  pressure  of  that  gas.  This  leads  to  the  following 
tissue  equations  for  and  that  are  used  in  the  model: 


dP, 


O. 


(28) 


(29) 


- - (P  -  P  - 

dt  0„  O,,  O,,  O,, 

2  2  2  2 


dP 


CO. 


- =+D  (P  -  P  )  +  M 

dt  CO^  CO^  CO^  CO^ 


■.-.'ich  being  i/RT,  the  ratio  of  nicle  concentration  to  partial  pre 

The  gas  transport  equations  24,  25,  28  and  29  are  solved  numerically 

with  the  blood  flow  equations  expressed  earlier  using  Hamming's  fourth 

order  predictor-corrector  integration  with  a  fourth  order  Runge-Kutta 
(37) 

starting  procedure.  This  formulation  consists  of  84  coupled,  nonlinear, 
simultaneous,  differential  equations  in  its  present  form  (including  the 
transport  of  only  and  CO^).  The  initial  results  are  given  in  the  next 
section  for  example  cases. 

RESULTS 


Although  the  model  has  been  tested  for  several  artificial  conditions 
(constant  blood  flow  with  constant  alveolar  gas  pressures,  constant  blood 
flow  with  sinusoidally  varying  alveolar  gas  pressures,  and  a  pumping  iieart 
with  constant  alveolar  gas  pressures),  results  will  be  presented  in  this  paper 
only  for  the  complete  model.  Specifically,  results  are  presented  here  for: 
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the  details  of  the  heart  pulse,  a  subject  exposed  to  a  sinusoidal  variation  of 
G^,  and  a  supine  subject  executing  passive  breathing. 

Figure  3  compares  the  details  of  physiologic  measurements  of  the 
heart  pulse  to  model  results.  The  model  was  run  for  6  seconds  of  real  time 
and  the  results  were  plotted  for  the  last  second  of  this  period.  Only  the 
artifacts  (0,C,0',C',a)  caused  by  valve  movements  are  absent  in  the  model 
results.  This  is  a  direct  consequence  of  the  omission  of  tissue  inertia  as 
discussed  previously.  As  it  is  not  our  purpose  to  model  every  detail  of 
heart  mechanics,  the  results  given  are  considered  to  be  in  satisfactory 
agreement  with  physiologic  measurements  for  proposed  applications  of  the 
model. 

The  differences  between  physiologic  values  and  model  values  for 
stroke  volume  are  due  to  a  heart  rate  of  70  beats/min.  in  the  model  while 
the  physiologic  measurements  were  taken  at  75  beats/min. 

Figure  4  shows  the  left  ventricular  pressure  on  an  expanded  scale  and 

includes  the  forcing  functions  as  specified  in  the  model.  Left  atrial 

pressure  has  also  been  plotted  so  that  the  effects  of  the  forcing  functions 

can  be  seen.  These  heart  pulses  were  generated  without  breathing  motion 

of  the  chest  wall  and  diaphragm  but  details  of  the  heart  pulse,  as  simulated 

by  the  model,  are  essentially  unchanged  by  passive  breathing  or  by  a 

sinusoidal  variation  of  G  of  amplitude  +  .5G  .  This  is  the  case  even 

z  -  z 

though  breathing  does  affect  blood  flows  through  changes  in  the  pressure 
within  the  chest  cavity,  and  hence  transmural  pressures  of  circulatory 
elements,  and  changing  G  does  affect  blood  pressure  throughout  the 
model  through  the  gravity  term  as  discussed  previously.  These  affects  are 
small  under  these  simulated  conditions,  but  will  become  more  significant  in 
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forced  breathing  and  for  larger  amplitudes  of  changes. 

Model  results  showing  capillary  blood  volumes  for  a  subject  exposed 
to  sinusoidal  variation  of  on  the  body  axis  with  amplitude  0.5g,  with 
positive  acceleration  being  directed  head-to-foot,  are  given  in  Figure  5. 
The  capillary  blood  volumes  resulting  from  the  model  run  without  this 
variation  of  G  are  shown  for  comparison.  Notice  that  the  rapid  changes  in 
capillary  volumes  in  the  upper  plots  are  small,  but  are  clearly  the  result  of 
the  ventricular  pressure  pulses.  A  phase  shift  can  be  noted  between  thr- 
heart  pulse  at  the  ventricle  and  the  volume  change  at  thp  ,.apillaries.  The 
symbol,  +,  in  the  figure  indicates  the  time  al  which  the  peak  of  systole  is 
reached.  A  similar  phase  shift  can  be  seen  in  the  variation  of  capillary 
volumes  due  to  the  applied  acceleration  in  the  lower  portion  of  Figure  5. 

Th^»  variations  in  capillaiy  volume  seen  here  are  due  to  both  the 

o 

variation  of  G^  and  the  ventricular  pulse.  The  180  phase  shift  between  the 
head  and  leg  capillary  volume  was  to  be  expected  as  gravitational  forces 
tend  to  reduce  (increase)  the  volume  of  blood  in  the  head  at  the  same  time 
they  tend  to  increase  (reduce)  the  volume  of  blood  in  the  legs. 

Partial  pressures  in  selected  chambers  of  the  model  are  shown  in 
Figures  6  and  7  for  a  supine  subject  executing  passive  breathing.  Although 
the  model  produces  flows,  pressures,  volumes  and  gas  concentrations  for  all 
chambers,  as  described  previously,  only  the  values  for  a  few  of  the 
chambers  have  been  selected  to  be  presented  here.  A  number  of 
observations  can  be  made. 

First,  the  partial  pressure  of  CO^  in  the  capillaries  is  almost  identical 
to  that  in  the  alveoli  whereas  the  partial  pressure  of  m  the  capillaries  is 
about  7  mm  Hg  below  that  in  the  alveoli.  This  results  because  the 
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diffusing  capacity  of  CO^  is  twenty  times  that  of  because  of  the 
different  molecular  weights  of  the  molecules. 

Secondly,  there  is  a  phase  difference  between  changes  in  partial 

pressures  in  alveoli  and  capillaries  for  both  and  CO^.  This  is  because  it 

takes  some  time  for  the  diffusion  processes  to  occur  across  the  vessel 

walls.  The  large,  periodic  drops  in  the  partial  pressure  of  in  the 

capillaries  shown  in  Figure  6  are  a  result  of  the  sudden  influx  of  veneous 

blood  resulting  from  ventricular  contractions.  This  would  be  reflected  in 

the  body  as  the  variation  in  the  distance  that  capillary  blood  must  travel 

through  the  lungs  (assuming  plug  flow)  before  it  becomes  fully  oxygenated. 

The  fluctuation  is  much  less  pronounced  in  as  shown  in  Figure  7  due 

CO2 

to  the  greater  diffusion  rate  of  CO^. 

The  fluctuations  in  the  partial  pressures  in  the  alveolar  blood 
chambers  are  damped  dramatically  as  the  blood  moves  from  one  chamber 
to  the  next.  This  effect  is  due  to  the  assumption  of  perfect  mixing  in  the 
model.  In  fact,  by  the  time  the  blood  reaches  the  arteries,  there  is  almost 
no  fluctuation  in  partial  pressures.  These  effects  are  more  pronounced  in 
the  model  than  the  real  system  because  perfect  mixing  does  not  occur  in 
the  real  system. 

A  slight  drift  is  noticeable  in  the  partial  pressure  of  and  CO^  in 
the  pulmonary  veins.  This  is  due  to  the  nonlinear  system  approaching  a 
slightly  different  equilibrium  under  pulsed  flow  as  compared  to  steady  flow. 
By  the  end  of  20  seconds  of  real  time,  however,  the  drift  of  values  has 
apparently  stopped  at  new,  dynamic  equilibrium  values.  These  values  are 
within  the  range  of  physiologic  measurements  as  shown  in  Table  II. 
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SUMMARY  AND  CONCLUSIONS. 

The  circulatory  model  described  here  appears  to  be  the  most 
complete  circulatory  model  published  to  date.  It  has  been  shown  to 
accurately  describe  pressures,  flows,  and  volume  changes  due  to  heart 
dynamics  in  response  to  a  simple  applied  forcing  function  simulating 
contraction  of  the  myocardial  muscles  (Figure  4  and  Appendix).  The 
transport  of  blood  throughout  the  model  subject,  produced  in  this  manner, 
has  been  shown  to  respond  to  small  variations  in  whole  body  acceleration  in 
a  realistic  manner.  In  addition,  the  circulatory  model  has  been  shown  to 
simulate  variations  of  gas  concentrations  at  the  mouth,  alveolar  chambers, 
and  throughout  the  circulation  when  interfaced  with  an  existing  lung 
model.  Gas  exchange  and  transport  has  been  modeled  in  such  a  way  that 
only  metabolic  rates  and  gas  concentrations  need  be  specified  to  provide 
predictions  of  temporal  variation  of  gas  concentrations  throughout  the  body 
under  passive  breathing.  The  detail  of  direct  physiological  correspondences 
included  in  this  cardiopulmonary  model  make  it  ideal  as  a  complete 
experimental  subject  for  physiological  studies.  Presently  the  model 
requires  about  25  seconds  of  CPU  time  on  a  Honeywell  6660  to  simulate 
one  second  of  real  time  but  this  can  be  reduced  on  a  faster  machine.  As 
was  mentioned  in  the  introduction,  only  the  details  of  the  circulatory 
portion  of  the  cardio-pulmonary  model  without  feedback  have  been 
presented  in  this  paper. 

The  circulatory  model  has  been  designed,  however,  so  that  feedback 
mechanisms,  external  body,  forces  arbitrary  breathing  maneuvers,  and 
breathing  of  different  gas  mixtures  can  easily  be  incorporated.  The  first 
stage  of  this  adaption  will  be  to  incorporate  neural  and  chemical  controls 
within  the  model.  This  is  an  essential  part  of  a  complete  circulatory  model 


27 


as  the  human  circulatory  system  responds  quite  differently  to  stimulus 
when  these  controls  are  blocked  as  compared  to  the  normal  response  with 
feedback.  Even  so,  it  is  important  that  a  circulatory  model  simulate 
supine,  resting  conditions,  even  without  feedback  mechanisms  in  operation, 
so  that  the  results  of  feedback  do  not  mask  deficiencies  in  the  primary 
model.  We  believe  that  this  has  been  shown  to  be  the  case  in  the  model 
presented  here. 

In  addition  to  the  incorporation  of  feedback,  the  model,  as  presented 
here,  must  be  modified  slightly  to  study  larger  whole  body  accelerations. 
The  modification  must  be  made  as  a  result  of  numerical  instabilities 
resulting  from  the  form  of  Equation  12,  which  models  the  pressure-volume 
relationship  in  the  circulatory  elements.  This  equation  predicts  an 
asymptotic  volume  at  infinite  pressurization  which  results  in  a  numerical 
instability  whenever  the  volume  of  a  circulatory  element  approaches  this 
asymptotic  value.  In  order  to  study  large  and  sustained  whole  body 
acceleration  where  some  vessels  will  fill  to  their  maximum  capacity, 
circulatory  element  volumes  in  these  vessels  which  approach  the  asympto¬ 
tic  volumes  must  be  handled  in  a  slightly  more  complex  numerical  manner 
than  that  presently  used.  Studies  to  incorporate  this  change  are  underway. 

Further  studies  with  the  model,  that  are  presently  being  considered, 

include:  sinusoidal  variation  of  larger  amplitude  and  a  variety  of 

(41) 

frequencies  for  comparison  to  experimental  data  of  Dr.  C.  F.  Knapp  ; 
linear  increases  of  g  to  large  sustained  values  for  normal  subjects  and 
subjects  with  borderline  abnormalities  of  the  cardiopulmonary  system;  and, 
ultimately,  simulation  ,  of  g  variation  encountered  in  aerial  combat 
maneuvers.  Currently  a  study  of  alteration  of  the  ventilation-perfusion 
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distribution  in  the  lungs  due  to  elevated  g  values  is  underway.  These  model 
studies  will  complement  experimental  studies  currently  in  progress  at  the 
U.  S.  Air  Force  School  of  Aerospace  Medicine.  Not  only  can  the  model  be 
used  for  these  studies,  but  it  also  can  be  used  as  a  tool  to  study  disease. 
For  example,  arterial  sclerosis,  emphysema  and  hypertension  are  only  a  few 
of  the  possibilities. 
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APPENDIX 

Forcing  Function  Used  in  the  Complete  Model  in  the  Heart  Chambers 

r-r-  -A 

FF  =  Be 

where  B  =  amplitude  of  the  pulse  and 

4  2  1 

A  =  l/2(t-t  -a/2)  /a  ,7>t-t  >0  where  o  =  pulse  rate  =  70/sec. 
0  <j  o 

at  t  -  t„  =  1/  6,  t  reset  to  t„=:  t. 

0  o  0 

0  =  (a/4)^ 
a  =  pulse  duration 


t^  =  starting  time  of  the  pulse 


u 

^0 

a 

(sec) 

(sec) 

(mml- 

Rt  Atrium 

.26 

.35 

5.0 

Rt  Ventricle 

.43 

.53 

20.0 

Lt  Atrium 

.26 

.35 

15.0 

Lt  Ventricle 

.43 

.53 

110.0 
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SECTION  B 

Pressure-Volume  Relationship 
in  Circulatory  Elements 


ON  THE  PRESSURE-VOLUME  RELATIONSHIP 


IN  CIRCULATORY  ELEMENTS* 

H.  H.  Hardy  anJ  R.  E.  Collins 
The  University  of  Houston,  Houston,  Texas  77004 

ABSTRACT 

An  empirical  equation  for  the  in-vivo  pressure-volume  relationship  in 
circul_tory  elements  is  proposed  as  a  universsil  form  valid  for  all  circulatory 
elements.  This  equation  contains  only  two  fundamental  physiologic  parameters, 
one  a  compliance  constant,  k,  and  the  other  a  maximum  element  volume,  V^.  It 
is  suggested  that  k  has  the  same  value  for  similar  elements  in  all  normal 
individuals  and  is  proportional  to  total  tissue  volume  of  the  subject.  Fits  of 
this  equation  to  published  atrial,  ventricular,  venous,  and  arterial  data  are 
presented  which  support  these  conjectures.  Combined  with  the  PoiseuUe 
relationship,  this  equation  is  shown  to  produce  a  fit  to  data  of  flow  resistance 
versus  transmural  pressure  for  the  pulmonary  vascular  system  of  a  dog. 
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Introduction 

Since  Bergel^^^  first  established  the  nonlinear  nature  of  vessel  elastance  a 
number  of  empirical  relationships  have  been  proposed  to  describe  the  dependence 
of  blood  volume  on  transmural  pressure  for  circulatory  elements.  The  relation¬ 
ships  introduced  by  Gessner^^^  Attinger^^^  and  Glantz^^^  provide  good  fits  to 
experimental  data  but  both  contain  a  multiplicity  of  parameters  which  vary  from 
vessel  to  vessel  with  no  apparent  physiological  basis  for  their  choice.  Gaasch^^^ 
has  proposed  the  simpler  form 

(1)  P  =  be’^'^ 

to  relate  transmural  pressure,  P,  to  end-diastolic  ventricular  volume,  V,  with  b 

and  k  empirical  constants.  This  relationship  requires  the  assignment  of  values  to 

(4) 

only  two  parameters  but,  as  pointed  out  by  Glantz  ,  it  is  subject  to  the  rather 
severe  criticism  of  predicting  zero  blood  volume  for  positive  transmural  pressure 
and  the  fit  to  data  deteriorates  as  the  accuracy  of  data  is  improved. 

In  this  brief  note  we  propose  a  new  form  of  equation  to  relate  transmural 
pressure  and  blood  volume,  valid  for  all  circulatory  elements  with  the  possible 
exception  of  capillary  beds,  which  is  free  of  these  criticisms. 

The  Proposed  Equation 

The  equation  we  propose  to  relate  transmural  pressure,  P,  and  blood 


1 


2 


volume,  V,  in  q  circulatory  element  can  be  deduced  by  the  following  argument: 
We  observe  that  data  of  V  plotted  versus  P  always  indicate  a  limiting  asymptote 
at  a  volume,  Vm,  characteristic  of  the  circulatory  element.  Thus  the  derivative 


(2) 


dV 

cfP 


=  f(V) 


can  be  viewed  as  a  function  of  V  which  approaches  zero  at  Vm  and  therefore  an 
expansion  of  f(V)  in  a  Taylor  series  in  (V  -  V^),  for  V  <  V^,  gives 


(3)  ^  =  k(V-V^^)  +  higher  order  terms 

with  k  the  value  of  -df/dV  at  V^.  Then,  neglecting  terms  of  higher  than  first 

order  in  (V-V  ) 
m 


(4) 


k  = 


1  dV 


m 

defines  k  as  a  constant,  having  dimensions  of  compliance,  characteristic  of  the 
circulatory  element.  Indeed,  k  characterizes  the  elastic  property  of  the  vessel 
wall  and  surrounding  tissue  while  the  other  parameter  in  this  equation,  V^, 
represents  the  limiting  volume  of  the  vessel  at  infinite  pressurization. 

We  note  that  since  the  right  member  of  Equation  (4)  can  be  multiplied  and 
divided  by  a  scale  factor,  say  a,  to  define  a  volume  V '  as  aV  having  a  maximum 
value  V  given  as  it  is  reasonable  to  expect  this  equation  to  apply  to 

circulatory  elements  of  different  sizes  with  one  value  of  k  applicable  to  all 
elements,  regardless  of  size,  if  these  are  of  the  same  type. 

Data,  taken  from  the  literature,  appear  to  support  these  conjectures. 
However,  before  we  examine  these  data  we  wish  to  stress  the  significance  of  the 
volume-scale  independence  already  noted  in  Equation  (4).  Since  all  norm*!! 
individuals  of  a  given  specie  have  sizes,  or  volumes,  of  the  various  elements  of 
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their  body  in  direct  proportion  to  their  total  mass  we  expect  that  the  value  of  V^, 

and  also  V  at  a  given  P,  should  be  proportional  to  the  size,  or  mass  of  the 

individual.  Indeed  such  "scaling"  of  anatomical  systems  and  parameters  has  been 

(6) 

observed  to  hold,  even  across  species'. 

The  pressure-volume  equation  implied  by  Equation  (4)  is  simply  obtained 
by  integration  as. 


(5) 


,  V  -V 

P  =  P  ~  X  V  ^-V~ 

°  ’'m  ^o 


with  V^,  being  the  values  of  V  and  P,  respectively,  at  an  arbitrary  point  on  the 
pressure-volume  curve.  We  wish  to  stress  this  fact  because  this  equation  actually 
contains  only  two  physiological  constants,  k  and  V^,  characteristic  of  the 
circulatory  elements.  We  now  show  how  well  this  equation  fits  in-vivo  data. 
Experimental  Validation 

(7) 

Figure  (1)  shows  data  from  Dahn  for  the  pressure-volume  relationship  in 
the  leg  arteries  of  men.  Figure  (2)  shows  these  same  data  normalized  on  the 


volume  scale  by  dividing  by  leg  volume  and  translated  vertically  on  the  V-V^  scale 
to  coincide  as  presented  by  Dahn.  This  translation  is  equivalent  to  changing  the 
reference  point  P^,  V^.  Also  shown  is  the  fit  of  Equation  (5)  to  these  data.  Note 
that  one  value  of  k  applies  to  all  three  individuals  and  blood  volumes  are 
proportional  to  total  tissue  volume. 

This  equation  was  also  fit  to  in-vivo  data  of  vessel  diameter,  D,  versus 

(8) 

transmural  pressure,  as  shown  in  Figure  (3)  using  data  on  dogs  from  Morris  ,  by 
2 

writing  V  as  ttD  L/4  and  assuming  the  length,  L,  a  fixed  constant.  Again  the  fit 
is  very  good  and  a  single  value  of  k  applies  to  both  animals.  The  scaling  of  blood 
volumes  by  animal  size  could  not  be  confirmed  here  because  animal  sizes  were  not 
reported. 
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I'iKiiro  1.  I’lolhysmogrnphic  daln  from  Dnhn'  showInR  thr  in-vivo  prcBSurc- 
volume  relationship  for  leg  arteries  of  men. 
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Figure  2.  The  curves  of  Figure  (1)  have  been  made  to  coincide  by  dividing 
element  volumes  by  tissue  segment  volume  and  then  translating  on 
the  vertical  scale.  The  fit  of  Equation  (5)  to  the  data  is  shown  by  the 
solid  curve.  The  parameter  k  has  the  same  value  for  all  three 
subjects. 


Figure  3. 
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A  fit  of  Equation  (5),  assuming  cylindrical  vessel  geometry,  is  shown 

(8) 

for  data  obtained  by  Morris  of  the  X-ray  determined  diameter- 


pressure  relationship  in  the  femoral  veins  of  dogs.  The  parameter  k 
has  the  same  value  for  both  dogs. 
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Fits  of  this  equation  to  atrial  pressure-volume  data  were  also  tested  by 
(9) 

using  data  from  Payne  as  shown  in  Figure  (4).  This  fit  to  diameter  versus 

3 

pressure  was  made  by  assuming  a  spherical  configuration  with  volume  it  D  /6. 
Clearly  the  fit  of  the  equation  to  data  is  better  than  the  reproducibility  of  the 
data.  Here  again  one  value  of  the  compliance  constant,  k,  applies  to  both 
subjects.  VVe  also  fitted  this  equation  to  some  in-vitro  data  of  ventricular  volume 

(5) 

versus  pressure  obtained  by  Glantz  as  shown  in  Figure  (5).  These  data  were  not 

very  reproducible  and  the  fit  of  the  equation  was  not  as  good  as  for  other  data, 

probably  because  this  was  not  in-vivo  data. 

As  a  final  example  showing  validation  of  the  proposed  pressure-volume 

relationship  we  exhibit  the  dependence  of  flow  resistance  on  transmural  pressure 

predicted  by  this  equation  and  compare  it  to  data.  Again  assuming  the  cylindrical 
2 

volume  TT  D  L/4,  with  L  fixed,  for  blood  volume,  we  use  the  Hagen-PoiseuUe 

4 

expression,  yL/2TrD  ,  for  flow  resistance,  R,  and  our  Equation  (5)  then  to  relate  R 
to  transmural  pressure  through  the  dependence  of  diameter  on  pressure.  The 
resulting  equation  can  be  written  as, 

(6)  R  =  ^ 

with  A  a  constant  depending  on  vessel  length  and  blood  viscosity  and  the 
dependence  of  V  on  P  given  by  Equation  (5).  If  Equation  (5)  is  substituted  here  the 
result  has  the  form 

-kP 

(7)  R  =  (a  -  be 

This  equation  was  fitted  to  data  of  flow  resistance  versus  transmural  pressure  for 
the  pulmonary  vascular  system  of  dog  lungs,  as  obtained  by  Borst^l®^  This  is  shown 
in  Figure  (6).  Clearly  the  fit  is  excellant  over  the  entire  range  of  pressure  even 
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Figure  4.  Fits  of  Equation  (5),  assuming  spherical  geometry,  to  in-vivo  end 

(9 

diastolic  atrial  pressure-diameter  data  in  dogs  obtained  by  Payne 
are  shown.  The  parameter  k  has  the  same  value  for  both  dogs. 


Figure  5.  Fits  of  Equation  (5)  to  in  vitro  ventricular  pressure-volume  data  of 
dogs,  obtained  by  Glantz^^\  are  shown.  The  parameter  k  has  the 
same  value  for  both  dogs. 


Figure  6. 
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fit  is  shown  of  Equation  (7)  to  data  obtained  by  Borst 


(10) 


of  flow 


transmural  pressure  in  the  pulmonary  vascular  system  in  a  dog. 


n 


though  the  pressure  flow  curve  in  the  lung  depends  on  several  physiologic  variables 
(e.g.,  lung  inflation,  alvelor  pressure,  vasomotor  tone,  etc.).  This  fit  provides 
further  confirmation  of  the  basic  validity  of  Equation  (5). 

Conclusion 

The  equation  proposed  here  to  describe  the  pressure-volume  relationship 
of  in-vivo  circulatory  elements  has  been  shown  to  fit  a  wide  variety  of 
experimental  data.  Although  a  direct  physiological  basis  has  not  yet  been 
established  it  appears  that  the  two  parameters,  k  and  V^,  of  the  equation  are^  in 
fact  physiological  constants  characteristic  of  circulatory  elements  and  that  the 
value  of  is  scaled  by  the  body  size  of  the  subject.  It  also  appears  that  the 
value  of  k  may  be  characteristic  of  each  anatomically  distinct  type  of 
c'rculatory  element.  Indeed  it  may  very  well  prove  possible  to  distinguish  an 
abnormal  circulatory  element  by  the  extent  to  which  the  k  value  differs  from  the 
universal  value  characteristic  of  the  element  type,  but  this  remains  to  be 
demonstrated. 
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ABSTRACT 


A  comprehensive  mathematical  model  of  pulmonary  mechanics 
and  regional  ventilation  is  developed  and  used  to  examine  the 
combined  effects  of  high  Whole  Body  Acceleration  (WBA)  and  altered 
body  position  on  respiratory  mechanics  and  ventilation.  The  model 
consists  of  a  system  of  driving  pistons  representing  the  respira¬ 
tory  muscles,  a  five  chambered  bronchial  tree,  and  six  parenchyma 
chambers.  It  also  simulates  gravitational  forces  on  the  masses  of 
the  chest,  diaphragm  and  abdomen  and  resulting  effects  on  fluid 
pressures  within  the  chest  and  abdominal  cavities.  Furthermore 
the  model  incorporates  the  effect  of  the  hydrostatic  pressure  gradient 
within  the  pleural  fluid  on  regional  lung  volumes  for  any  gravitational 
field.  The  model  is  shown  to  be  adequate  for  simulation  of  many 
clinical  pulmonary  function  tests.  Several  studies  are  conducted  to 
determine  effects  of  the  high  WBA  environment  on  respiration  and  how 
altering  body  orientation  relative  to  the  G  field  may  improve  G 
tolerance.  As  a  result  of  these  studies,  it  is  concluded  that  the 
optimum  seatback  angle  for  maintaining  near  normal  ventilation  under 
WBA  is  between  60°  and  75°. 
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I.  INTRODUCTION 

The  effects  of  artificially  altered  gravitational  forces 
caused  by  whole  body  accelerations  (WBA)  on  the  human  respiratory 
system  were  first  recognized  as  a  problem  for  study  shortly  after 
the  development  of  modern  aviation.  Today,  with  the  advent  of  high 
performance  aircraft  in  use  by  our  military,  as  well  as  similar  planes 
in  use  by  other  countries  around  the  world,  this  problem  has  become 
a  critical  limiting  factor  in  the  training  of  pilots.  While  some 
planes  in  use  today,  specifically  the  F-15,  YF-16,  and  YF-17,  have  the 
capability  of  withstanding  up  to  16-G  maneuvers,  a  pilot  will  usually 
blackout  before  12-G's.  As  has  been  demonstrated  by  Burtonf"'^  even 
a  small  advantage  in  WBA  tolerance  limits  will  almost  certainly  give 
a  pilot  the  needed  superiority  during  an  aerial  combat  engagement. 
Consequently,  any  means  of  defining  the  physiological  parameters 
involved  in  pilot  failure  and,  hopefully,  how  these  might  be  modified 
to  improve  pilot  tolerance  of  WBA,  prove  to  be  of  scientific  interest 
and  in  the  interest  of  our  national  security. 

Thus  far,  the  primary  tool  for  investigating  the  effects  of  WBA 

has  been  the  human  centrifuge.  The  history  of  this  device  dates 

(2) 

back  at  least  to  1815'  '  when  a  centrifuge  was  used  for  the  treatment 

(3  ] 

of  mental  patients;  however,  as  pointed  out  by  Gauer'  '  in  his  review 
of  the  history  of  the  machine,  it  had  its  maximum  healing  effect  only 
when  a  second  treatment  was  suggested.  The  centrifuge  consists  of 
a  gondola  on  a  long  arm  which  may  be  spun  at  controlled  speeds.  (See 
Figure  1).  WBA  is  given  as  multiples  of  earth  normal  gravity  (G's) 


and  is  calculated  as 
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(1) 


WBA  = 


2/  ^2 
4Tr  (ros)  r 


where  rps  is  revolutions  per  second,  r  is  radius,  and  g  is 
the  gravitational  constant.  The  nomenclature  used  to  define  the 
direction  of  forces  generated  are  given  in  terms  of  body  geometry 
(See  Figure  2).  The  centrifugal  field  must  be  summed  vectorially 
with  Earth's  gravitational  field  to  obtain  the  resultant  field.  To 
eliminate  confusion,  all  future  references  to  WBA  will  be  assumed 
to  refer  to  the  magnitude  of  the  resultant  force  and  direction 
indicated  by  degrees  off  the  axis,  i.e.,  is  0°,  is  90°, 

-G^  is  180°,  and  -G^  is  270°.  The  axis  along  the  resultant  force 
line  will  be  referred  to  as  the  G-axis.  Left  and  right  lateral 
forces  will  not  be  examined  at  this  time. 

While  the  centrifuge  remains  the  primary  tool  used  to  investigate 
the  effects  of  WBA,  it  has  many  drawbacks  and  data  collected  on  this 
device  in  the  study  of  the  pulmonary  system  is  sketchy  at  best.  Many 
experiments  cannot  be  run  on  animals  as  they  require  cooperation  of 
the  subjects.  When  working  with  human  subjects,  for  the  most  part, 
only  noninvasive  techniques  are  allowed.  Experimental  control  is 
difficult  as  fatigue  sets  in  rapidly,  limiting  the  data  collected 
on  a  subject  during  a  single  session.  Work  with  human  subjects  is 
further  complicated  by  the  fact  that  high  WBA  studies  are  always 
associated  with  a  certain  amount  of  risk  and,  besides  the  parameters 
of  interest  to  the  study,  many  vital  signs  must  be  monitored  for 
which  medical  personnel  must  be  present.  It  is,  therefore,  desirable 
that  a  good  model  of  the  pulmonary  system  be  developed  which  may  be 
used  to  direct  and  augment  studies  using  the  human  centrifuge. 
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Many  models  of  the  pulmonary  system  have  been  developed  in 

recent  years.  One  of  the  most  frequently  encountered  is  the  single 

chambered  model.  A  very  detailed  model  of  this  type  was  developed 

f4) 

by  Suwa  and  Bendixin^  '  which  incorporated  diffusion  and  perfusion. 

(5 ) 

Scheid  and  Piiper'  '  employed  such  a  model  to  evaluate  the  effect 
of  respiratory  patterns  of  gas  washout  data.  Among  those  to  apply 
such  a  model  to  the  study  of  lunq  mechanics  are  Fry  and  Hyatt^^^ 
who  examined  the  effects  of  tissue  elastance  and  pleural  pressure 
on  flow  resistance  in  a  single  airpassage  and  Collins,  et  al.,^^^ 
whose  model  differed  from  that  of  Fry  and  Hyatt  in  that  it  incorporated 
the  compliance  of  the  respiratory  muscles  such  that  the  role  of  the 
active  components  of  the  ventilation  system  could  be  examined.  The 
last  two  models  are  among  those  most  firmly  based  on  sound  physical 
principles.  Features  from  each  are  incorporated  into  our  model  as 
discussed  in  the  next  section. 

A  few  multicompartment  models  have  been  advanced  for  the 
evaluation  of  ventilation  tests.  Gomez,  et  al .  ,  ^ ^  have 

developed  a  very  sophisticated  mathematical  treatment  of  such  a 
model,  but  it  has  been  criticized  for  assuming  all  chambers  to  be 
ventilated  in  paral lei . Both  Saidel,  et  and 

West^^^’^^’^^^  have  developed  similar  models  which  included  both 
series  and  parallel  ventilation,  which  is  a  more  realistic  approach, 
and  is  incorporated  into  our  model.  Their  models  were  designed  to 
study  ventilation,  however,  and  do  not  lend  themselves  readily  to 
the  study  of  mechanics. 
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Among  those  models  designed  specifically  for  the  study  of 
mechanics  are  those  of  Head,  et  al.,^^^^  and  Lambert  and  Wilson, 
both  of  which  were  designed  to  study  elastic  characteristics  of  the 
lung.  Their  work  may  be  criticized  as  the  physical  laws  employed 
are  valid  only  for  small  deformations,  thus  entailing  linear  theory. 

Two  other  types  of  model  commonly  encountered  in  the  literature 
are  (a)  those  which  employ  electrical  analogs,  noteworthy  of  which 
are  those  of  Grimby,  et  al.[^*^^  and  Hilberman,  et  al.l^^^  and 
(b)  those  which  apply  curve  fittina  techniques  to  clinical  and 
experimental  data  and  attempt  to  correlate  results  with  physiologic 
parameters.  The  first  of  these,  by  virtue  of  its  inherent  linearity, 
can  only  be  addressed  to  linear  portions  of  lung  dynamics  which 
restricts  its  usefulness  to  a  narrow  range  around  passive  breathing. 
Models  of  the  second  type  are  among  those  most  commonly  encountered 
in  the  literature;  however,  they  yield  virtually  no  direct 
information  on  lung  mechanics. 

One  model  of  note,  which  constitutes  a  separate  class,  is 

(22) 

that  advanced  by  Kilpper  in  ris  doctoral  dissertation;  unfortunately, 
it  has  never  appeared  in  public  literature.  His  model  consists  of 
a  con'*'inuous  lung  tissue  system  rather  than  a  chambered  system. 

While  this  model  is  useful  in  the  study  of  regional  ventilation,  it 
does  not  incorporate  a  realistic  model  of  the  surrounding  tissue. 

’’hough  not  complete,  this  list  is  representative  of  what  has 
been  accomplished  in  this  field.  All  of  these  models  suffer  from 
one  or  more  of  the  following  restrictions  to  their  usefulness; 
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i)  They  rely  too  heavily  on  model  parameters  which  are  not 
clearly  identifiable  with  anatomical  or  physiological 
elements . 

ii)  They  assume  linearity  of  such  parameters  as  airway 

resistance  or  compliance  when  nonlinearity  ■'s  known  to 
exist. 

iii)  They  assume  that  variables  such  as  gas  concentrations  or 
alveolar  pressures  are  uniform  throughout  the  system. 

( 23 1 

As  was  pointed  out  by  Hyatt  and  Black^  '  in  their  recent  review 
article,  it  will  not  be  until  an  adequate  mathematical  model  of  the 
pulmonary  system  is  developed  that  the  various  pulmonary  function 
tests  now  available  can  be  ful ly  exploited  as  diagnostic  tools.  They 
qualify  an  "adequate"  model  as  a  mul tichambered  model  based  upon 
sound  physical  principles  -  rather  than  curve  fitting  -  which  incorporates 
the  basic  nonlinear  nature  of  the  individual  p'ji'rionary  elements.  These 
comments  can  be  extended  to  say  that  if  a  model  is  not  caoable  of 
simulating  the  various  pulmonary  function  tests,  it  most  certainly 
will  not  be  of  use  in  examining  the  effects  of  WBA. 

This  dissertation  will  describe  the  continued  development  of 

a  comprehensive  numerical  simulator  of  the  human  pulmonary  system 

first  advanced  by  the  author  as  partial  fulfillment  of  the  requirements 

(24 ) 

for  the  masters  degree'  and  its  application  to  specific  effects 
of  WbA. 

Chapter  II  will  provide  a  description  of  the  analomical 
approximations  and  mathematical  relations  developed  to  portray  the 
system.  This  will  be  presented  much  as  it  was  in  the  above  mentioned 


masters  thesis. 
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Chapters  III  and  IV  will  describe  recent  modifications  that 
have  been  made  to  allow  the  model  to  depict  adequately  'rway 
resistance  and  lunq  tissue  compliance  respectively.  Soi:,e  preliminary 
results  indicating  improved  correspondence  between  model  predictions 
and  observed  data  are  presented  in  these  chapters,  the  most  noteworthy 
of  which  are  the  pulmonary  pressure-volume  curve  and  the  graded 
effort  maximum  expiratory  flow  volume  curves. 

Chapter  V  will  examine  the  effect  of  VJBA  on  the  respiratory 
system  through  four  studies  utilizing  the  pulmonary  simulator,  and 
Chapter  VI  will  relate  these  studies  to  existing  knowledge  in  such 
a  way  that  an  optimum  seatback  angle  for  relieving  the  adverse  effects 
of  WBA  may  be  defined. 
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II.  THE  MODEL 

A  matheinatical  model  of  the  human  pulmonary  system  has  been 

developed  which  incorporates  the  desirable  features  of  many  models 

previously  developed  and  satisfies  the  requirements  for  a  model  set 

forth  by  Hyatt  and  Black  (see  page  7  ).  It  incorporates  the  volume 

dependent  nonlinear  resistance  described  ty  Fry  and  Hyatt^^^,  the 

'15  16  17) 

series  and  parallel  ventilation  concept  advanced  by  West'  ’  '  '' 

and  Saidel,  et  and  the  compliance  of  respiratory  muscles 

included  in  Coll'ins'  model^^^  It  also  includes  several  parameters 
such  as  the  spatial  distribution  of  mass  of  surrounding  deformable 
tissue  and  the  hydrostatic  pressure  gradient  in  the  pleural  fluid 
which  have  been  overlooked  in  previous  models,  though  essential  to 
the  study  of  WBA. 

The  model  advanced  by  Collins,  et  al.l^^  consists  of  a  single 
driving  piston,  representing  the  muscles  of  respiration,  which  acts 
across  a  fluid  filled  pleural  space  on  an  alveolar  chamber  which  is 
connected  to  the  outside  of  the  body  via  a  single  elastic  passage 
of  the  type  employed  by  Fry  and  Hyattl^^  Though  adequate  for 
examining  gross  lung  measurements,  such  as  pulmonary  function  tests 
under  specific  conditions,  this  model  can  not  be  used  to  explain 
results  in  terms  of  regional  characteristics.  Further,  sin^e  all 
elastic  tissues  exterior  to  the  lungs  were  lumped  into  a  single  piston, 
irrespective  of  geometry,  their  model  proved  inadequate  in  predicting 
results  under  chanqing  conditions  such  as  whole  body  acceleration  or 
charging  body  orientation  with  respect  to  the  G-axis. 
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The  respiratory  muscles  are  scattered  throughout  the  trunk 
and  neck  regions  of  the  body  ana  do  not  act,  nor  are  they  acted 
upon,  by  gravity,  as  a  single  unit.  They  should  be  subdivided  into 
a  minimum  of  three  anatomical  groups,  namely: 

(1)  thoracic  muscles  -  muscles  which  act  on  the  rib  cage 
and  spi . e , 

(2)  abdominal  muscles,  and 

(3)  diaphragm  muscles. 

Lung  volume  is  affected  by  the  weight  of  each  of  these  groups 
according  to  the  architectural  design  of  the  human  body  and  its 
orientation  in  the  gravitational  field.  Specifically,  lung 
volume  would  be  most  affected  by  the  mass  of  the  (a)  diaphragm 
when  the  G-axis  extends  through  the  cephalocaudal  axis,  6^  axis, 

(b)  the  abdominal  muscles  when  the  G-axis  passes  through  the 
ventrodorsal  axis,  axis,  and  (c)  the  chest  when  the  G-axis  is 
perpendicular  to  the  ribs  so  as  to  maximize  the  moment  arm  of  force 
on  the  ribs.  Since  the  diaphragm  and  abdomen  act  in  series  and 
since  the  diaphragm  acts  only  during  inspiration  while  the  abdomen 
acts  only  during  expiration,  these  may  be  thought  of  as  a  single 
group  if  the  necessary  geometry  is  included  to  describe  the 
above  mentioned  mass  distribution  effects.  This  allows  the  two 
piston  model  depicted  in  Figure  (3)  where  the  geometry  of  mass 
distribution  is  illustrated  in  Figure  (4). 

The  lungs  are  suspended  in  the  thoracic  chamber  with  points  of 
connection  with  the  rest  of  the  body  being  only  the  hilus  and 
Dulmonary  ligaments,  through  which  the  trachea  and  pulmonary  vessels 


abdominal  cavity 
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pass.  It  is  thus  justifiable  to  treat  the  lunqs  as  a  separate 

organ  though  they  cannot  function  as  a  separate  unit  as  they 

(25 ) 

contain  no  internal  driving  mechanism,  i.e.,  muscle  tissue;  '  The 
ventilation  portion  of  this  organ  consists  of  (a)  an  air  delivery 
system,  the  conductive  zone,  and  (b)  the  functional  or  respiratory 
zone. 

{  ?fi  ) 

The  conductive  zone  was  extensively  examined  by  Wiebel,  who 
developed  two  alternate  morphological  models  of  the  "respiratory 
tree".  He  meticulously  traced  and  measured  branches  of  this  system 
in  five  normal  human  lungs  for  the  first  five  branches  and 
statistically  for  all  remaining  branches  out  to  their  terminus.  His 
model  A  consists  of  the  average  values  for  all  measurements  for 
individual  segments  of  air  passages  and  is  characterized  as  an 
irregularly,  dichotomously  branching  system  of  airways  that  is 
continuous  for  between  sixteen  and  twenty-five  generations  of 
branches.  Model  B  is  a  regularly  branching  air  passage  system  in 
which  the  average  measurements  for  all  air  passage  segments  of  the 
same  generation  are  used.  Because  of  its  symmetry,  model  B  is  most 
often  used  by  mathematical  modelers.  The  model  herein  developed 
does  not  require  symmetry  and  either  of  Weibel's  models  may  be 
represented . 

After  an  average  of  nineteen  generations,  designated  as  the 
trachea,  bronchi,  bronchioles,  and  terminal  bronchioles  respectively 
according  to  their  relative  depth  within  the  system,  the  conductive 
air  passages  empty  into  the  acini  which  consist  of  the  respiratory 
bronchioles,  alveolar  ducts,  atrium,  and  alveolar  sacks,  all  of  which 
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are  lined  with  tiny  alveoli.  There  are  within  the  lungs  approximately 
300  x  10^  alveoli^^^^  and  it  is  only  through  the  walls  of  the  alveoli 
that  gas  exchange  between  air  and  blood  takes  place.  The  air  volume 
close  enough  to  these  specialized  structures  for  mixing  and  diffusion 
to  bring  it  into  contact  with  their  walls  has  been  designated  the 
parenchyma,  the  functional  unit  of  which  is  the  acinus. 

It  is  the  parenchyma  which  has  been  designated  the  alveolar  volume 
of  all  compartmental  models  including  the  one  herein  advanced.  The 
difference  between  this  volume  and  the  total  lung  volume  is  the 
anatomical  dead  space'  ^  and  corresponds  to  the  total  air  passage 
volume  within  the  model. 

The  parenchyma  may  -  and  has  been  -  subdivided  into  any  number 

of  individual  functional  units  for  modeling  purposes.  In  the  ideal 

model  this  number  would  correspond  to  the  number  of  acini  in  the  real 

lung  which  would  be,  since  they  occur  after  approximately  nineteen 
1 9  5 

branches,  2  or  about  5x10  ,  which  is  clearly  to  large  a  number 
for  modeling  purposes.  The  field  of  engineering  (a  discipline  that 
uses  models  extensively  and  from  which  many  of  the  models  of  biological 
systems  have  come)  emphasizes  simplicity  in  modeling  and  it  is  upon 
this  rule  that  the  single  chambered  models  are  based. 

An  effect  which  cannot  be  examined  with  a  single  chambered  model  and 
which  has  not  been  included  in  any  mul tichainbered  model,  to  the  knowledge 
of  this  investigator,  is  that  of  the  pleural  pressure  gradient  in  the 
thoracic  cavity.  Given  the  pleural  pressure  at  the  top  of  the  lungs 
with  respect  to  the  G-axis,  P  ,  ,  the  pleural  pressure  at  any 

P  0 

depth,  h,  along  the  G-axis  may  be  expressed  as 
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Ppl=Ppl^^qgh  , 

where  q  is  the  mass  density  of  lung  tissue  and  g  is  the  gravitational 
acceleration. 

The  effect  of  this  pressure  gradient  on  the  lungs  was  clearly 

(31 1 

demonstrated  by  Glazier,  et  al.,  '  in  his  examination  of  greyhound 

lungs  that  were  fixed  in  situ.  He  found  that  apical  and  basal 

alveoli  were  the  same  size  when  fixed  in  the  same  horizontal  plane, 

while  apical  alveoli  were  four  times  larger  than  basal  alveoli  when 

the  lungs  were  fixed  vertically.  This  gradient  in  unit  volume  has 

been  attributed  to  the  geometric  mass  distribution  around  the  lungs 

(32  33  341 

and  the  hydrostatic  pressure  gradient  in  the  pleural  fluid. ^  ' 

The  parenchyma  has  been  subdivided  within  the  model  (sec 
Figures  (3)  and  (5))  into  six  chambers  arranged  three  to  a  side  and 
spaced  at  distances  h.| ,  ^2^  and  h^  along  the  cephalocaudal  axis  and 
h^  and  h^  along  the  right  to  left  axis  in  the  coronal  plane.  The 
symmetry  imposed  by  this  geometry  tends  to  simplify  the  examination 
of  the  effects  of  the  hydrostatic  gradient  mentioned  above  by 
subdividing  each  side  of  the  lung  into  hydrostatic  levels,  thus 
allowing  the  simulation  of  any  body  orientation  with  respect  to  the 
G-axis . 

It  should  be  noted  that  this  breakdown  is  purely  artificial  and 
bears  no  relation  to  the  five  lobular  sections  of  the  lungs.  Strict 
adherence  to  anatomical  structure  may  be  imposed  by  effectively 
eliminating  one  alveolar  segment  on  the  left  by  setting  its  volume 
to  zero  and  the  resistance  into  it  to  infinity  and  adjusting  the  size 
and  relative  location  of  the  other  chambers  to  a  geometry  matching 
the  true  anatomical  state. 


The  airpassages  have  been  subdivided  into  a  system  of  five 
bifurcating  air  chambers  allowing  both  series  and  parallel 
ventilation  of  the  six  alveolar  chambers  as  indicated  in  Figures 
(3)  and  (5).  Within  the  framework  of  Weibel's  models,  Chamber  A-j 
in  Figure  (5)  corresponds  to  the  nasal  or  oral  cavities,  larnyx, 
trachea,  and  primary  bronchi.  Chambers  A2  and  A^  correspond  to  all 
secondary  and  possibly  tertiary  bronchi  in  the  lungs  and  all  lower 
order  bronchi  and  bronchioles  supplying  the  upper  third  of  the 
lungs.  Chambers  A^  and  A^  would  then  correspond  to  all  lower  order 
bronchi  and  bronchioles  supplying  the  lower  two-thirds  of  the  lungs. 

Measurements  of  in  vivo  human  lung  volume  cann^^t  be  made 
internal  to  the  lungs  but  are  generally  accomplished  indirectly  by 
a  family  of  instruments  called  spirometers While  these  vary 
greatly  in  design  they  all  serve  the  same  basic  purpose  -  to  measure 
the  volume  of  air  inspired  and  expired  by  the  subject  during  various 
breathing  maneuvers.  Chamber  Spr  in  Figure  (4),  which  completes 
the  gas  volume  portion  of  the  model,  represents  just  such  a  device. 
When  the  initial  volume  of  this  chamber  is  set  to  a  relatively 
small  value  (but  larger  than  the  volume  that  can  be  inspired  by  the 
lungs)  it  represents  an  ideal  spirometer,  meaning  that  it  is 
frictionless  and  operates  at  constant  pressure  and  temperature.  If 
the  volume  of  this  chamber  is  set  to  a  very  large  value  it  simply 
represents  the  atmosphere.  Source-sink  relations  for  individual 
gasses  are  included  within  this  chamber  such  that  it  can  be  used  to 
model  specialized  life  support  systems  such  as  are  found  in  a  pilot's 
cockpit  or  a  diving  sohere. 
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One  assumption  has  been  made  with  regard  to  the  air  chambers 
to  simplify  the  mathematics  involved  which  requires  some 
clarification;  all  air  chambers  in  the  model  are  considered  "ideal 
mixina'  chambers.  This  concept  will  be  explained  for  the  airpassages 
and  the  explanation  may  be  extended  to  the  other  chambers.  When 
air  (or  fluid)  passes  through  a  tube  it  travels  not  as  a  square 
wave  (plug  flow)  pushinq  all  old  air  in  front  of  the  new,  but  as  a 
paraboloidal  front  where  velocity  is  greatest  at  the  center  and 
may  be  zero  near  the  walls.  This  means  that  the  composition  of 
air  near  the  wall  will  closely  resemble  that  of  the  old  resident 
gas  while  air  near  the  center  will  reflect  that  of  the  inflowing 
gas.  This  analogy  becomes  somewhat  muddled  around  points  of  high 
resistance  where  turbulence  may  occur,  but  the  message  is  the  same  - 
composition  of  air  within  a  chamber  must  be  represented  as  some 
average  of  resident  and  inflowing  gas.  This  concept  has  been  well 
developed  by  Pedley,  et  al.'  '  The  assumption  herein  applied  is 
that  incoming  gas  instantly  mixes  with  all  resident  gas.  Considering 
that  instantly  here  means  within  one  finite  integration  step  (0.01  - 
0.001  sec),  that  the  volumes  of  the  chambers  are  only  0.05  -  0.20  L, 
and  that  the  ^low  rate  may  exceed  1100  L/min.,  this  is  a  safe 
assumption. 

The  development  of  the  equations  describing  this  system  starts 
with  the  assumption  that  the  volume  of  the  thoracic  cavity  is  equal 
to  the  air  volume  plus  a  constant,  K,  which  is  the  volume  of  pleural 
fluid  and  the  volume  of  lung  tissue,  i.e.. 
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(3)  = 

where  V-j.^  is  the  volume  contribution  of  the  thorax  ’  ^DA  the 
volume  contribution  of  the  diaphraqm-abdomen  comolex,  and 
for  j=l  to  6  are  the  six  alveolar  chamber  volumes. 

The  force  balance  equations  describing  the  pressures  acting 
across  the  muscle  and  parenchyma  chamber  walls  (see  Figure  6)  were 
developed  from  Newton's  second  law.  Assuming  the  inertial  terms 
to  be  negligible,  the  force  balance  equations  are: 


(4) 


Tx  ■  ^bs  ■  "^pl 


+  P.  +  P.1 

g  el 


Tx 


for  the  thorax, 


(5) 


DA 


^^bs  =  Ppl 


^b  % 


+  P 


el 


DA 


for  the  diaphragm-abdomen  complex,  and 


(6) 


pi 


+  P  +  P.1 
9  •  el 


for  each  of  the  six  alveolar  chambers,  j  =  1  to  6. 

Equations  describing  the  volume  dependence  of  the  elastic 
terms,  P^i ,  will  be  discussed  in  a  later  chapter.  The  basic  equatic 
of  pulmonary  statics  as  described  by  the  model  are  complete  when 
the  force  balance  equations  (Equations  4-6)  and  the  conservation 
of  pulmonary  volume  equation  (Equation  3)  are  coupled  with  the 
equation  of  state  for  the  gas  in  each  chamber,  i.e.,  assuming 
isothermal  conditions, 
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where  n  is  the  total  moles  of  gas,  R  is  the  ideal  gas  constant, 
and  T  is  temoerature. 

Two  additional  equations  are  necessary  to  describe  pulmonary 
dynamics;  a  mass  balance  equation  describing  the  gas  in  each  chamber 
and  a  flow  equation  describing  the  transfer  of  gases  between 
chambers.  The  latter  of  these  will  be  treated  in  a  later  chapter 
on  airway  resistance.  If  the  initial  concentration  of  gas  1  in 
chamber  j  is  known,  the  moles  of  gas  present,  n.jj,  may  be  determined 
by 


(8) 


^Ij  "  ^0 


Ij 


dn 


Ij 


dt 


dt 


v/here  n^,  is  the  initial  value  of  n,  ■  and 
^Ij 


(9) 


dn,  . 


Where  i  is  any  chamber  having  flow  into  Chamber  j,  k  is  any 
chamber  into  which  flow  in  passing  from  j,  T.^  is  the  volume  flow 
rate  from  i  to  j ,  and  C-j  ■  is  the  molar  concentration  of  gas  1 
in  chamber  i,  Diffj  is  a  function  reorescnting  gas  exchange  with 
the  blood  and  would  only  appear  in  the  alveolar  equations. 

The  complete  set  of  equations  consisting  of  the  derivative  of 
the  volume  equation  (Equation  3),  the  force  balance  equations 
(Equations  4-6),  the  equation  of  state  (Equation  7),  and  the  mass 
balance  equations  (Equation  8)  for  each  chamber  and  each  gas  as 
required,  constitutes  a  set  of  150  simultaneous,  nonlinear 
aifferential  equations.  (For  derivation  see  Appendix  A.)  These 
are  listed  in  Tab‘'es  I,  II,  and  III.  The  only  equations  listed 


//A  t/EOLAR 


EQ.UAT/ONS 


Ez  o.  /  Eo. 


E:  Eion 


/  -  ^ 


^,- 


EiE^L 


cy ' 


_ _. 

><  ^ 
i4- 


i/t-,  ~  .'^  ■'"'-X 

V 


7 


O' 


OA 


H  >;  -  ^  ''■,  '  R 


c»  ‘'/2A 


"  ^.,; 


■  O 


l5 


c^-J/ 


/n 

'^Jp  /r? 


/  4 


V 


r?uz 


7 


,  <^CoJ.'  — 

s:?-'  7 


7  7 

y  ^ 


33-6/ 


C, 


V 


/  /  ‘  • 

'~/  (  -  Cn  ■  I/: 

i<-  ^  -V  ■</  V  • 


(^>2  67 


6:^'. 


C./J-  V 


68-73 


3 


rY^j  ’■7,. 


V  7  ^ 


/ 


7V-  79 


^‘‘Y  ^<^/n  8^mE 


Equations  are  numbered  cumulatively  based  upon  the  chamber  number 
j»  7-12  and  gas  number  !■  1-4. 


*Set  Appendix  A  for  definition  of  m.  Roughly,  It  is  the  chamber  from 
which  gas  is  flowing. 


TRACHEO-  BRONCHIAL  EQUATIONS 


hcj.  No.  ^ 

^ Cf  us'i  /  on 

/-ZO 

4-  - 

2/ -  ^0 

% 

V/  - 

P-  = 

ZTA- 

£ 

J 

( 

C.J.  )  K- 
/•pr  ^E=-  N  ^ 

V5- 

- 

r/  -  rT 

% 

S' 

fr6  -66> 

% 

S'  V  f  -  S 

Equations  are  numbered  cumulatively  based  upon  the  chamber 
J*“  2-5  and  gas  number  1*  1-4 

table  JI 


5PfR0/^E  TER  Ea  UA  T /ONS 


Fa.  No. 


9 


udi:  lor^ 


/- V 


5" 


^  -- 


'Jl/ 


/?-r 


—  Cuj/ 


6 


^CJ/  '  C<.J,  R/ 


7-/0 


^jg/  ~  1/  C  '^jif  ^j/  ^ 

/ 


// 


S^J/  -  Cujrtf  N 


*4.7/  ^  UJ  ^ 


/ABLE  HI 


25 


which  have  not  been  discussed  thus  far  are  those  for  the  molar 
water  concentration,  ,  the  moles  of  water  in  a  chamber,  n^  , 

j  j 

and  the  source  term  for  water  vapor,  S  .  Since  water  vapor, 
unlike  the  other  gases  present,  has  a  saturation  point  dependent 
on  the  P  and  T  of  the  chamber,  a  source  term  must  be  supplied 
in  the  mass  balance  equation  for  water  to  allow  for  evaporation 
and  condensation.  For  simplicity  all  chambers  in  the  model  are 
assumed  to  be  at  the  saturation  point  at  all  times. 

The  equations  in  these  three  tables  constitute  the  complete 
set  of  differential  equations  required  to  describe  the  dynamics  of 
this  model.  These  were  algebraically  rearranged  such  that  they 
could  be  sequentially  solved  (see  Appendix  A)  and  were  programmed 
for  numerical  integration  on  a  high  speed  Computer  in  Fortran  IV 
language. 

The  model  is  always  initialized  in  the  relaxed  state  where 

and  Pp^  are  both  zero.  In  this  state  the  system  will  take 
on  its  equilibrium  volume,  which  is  the  functional  residual 
capacity  (FRC),  and  all  pressures  within  the  closed  air  system  will 
equal  p5p^-  P5pf.  f’bg  i^ust  be  specified  constants  or  time 

dependent  functions  representing  the  environmental  conditions  being 
modeled.  Given  the  initial  gravitational  acceleration  and  its 
orientation  with  respect  to  the  body  coronal  plane,  all  hydrostatic 
pressures  in  the  pleural  '^nace  may  be  calculated  using  Equation  (2) 
and  pressures  due  to  the  weight  of  each  muscle  group  may  be  calculated 
from  input  parameters  related  to  the  mass  of  me  structure  and  the 
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area  over  which  it  operates.  P^g  may  be  calculated  from  an 
equation  similar  to  Equation  0.  Since  the  elastic  pressures  are 
modeled  as  functions  only  of  volume,  Equations  (3-6)consti tute  a 
set  of  nine  equations  in  nine  unknowns,  namely:  Vj^,  Pp^ , 

and  Vj  for  j  =  1-6. 

To  solve  this  set  of  equations  for  initial  values  of  the 
dependent  variables  listed  above  a  Newton-Raphson  iteration  was 
first  tried.  The  equations  proved  to  be  very  oscillatory,  however, 
and  convergence  could  not  be  achieved.  The  modified  Newton-Rapson 
iteration  procedure  advanced  by  R.  P.  Donachi,^  '  in  which 
original  estimates  are  corrected  by  only  half  the  calculated 
correction  values,  was  attempted  but  it  was  found  that  either 
convergence  was  too  slow  or  the  system  would  converge  to  a  constant 
error.  We  further  modified  Donachi's  method  by  using  the  full 
correction  value  if  the  calculated  error  for  any  variable  did  not 
change  sign  during  three  consecutive  iterations,  returning  to  the 
half  step  any  time  sign  change  was  noted.  This  allowed  convergence 
within  an  acceptable  number  of  iterations. 

A  fourth  order  Hamming  predictor-corrector  algorithm  using  a 

(37 1 

fourth  order  Runge-Kutta  starting  p'-ocedure  '  was  chosen  as 
integration  technique;  however,  upon  attempting  to  run  the  mouel 
instabilities  were  indicated.  After  performing  a  stability 
analysis^  ’  it  was  found  that  the  oressure  equations  in  the 
airways  were  absolutely  unstable  and  the  problem  could  not  be 
corrected  using  the  equations  as  formulated. 
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To  correct  this  problem  we  assumed  a  steady  state  distribution 
of  pressures  in  the  airways;  that  is,  the  pressures  within  the 
airway  chambers  are  assumed  to  readjust  instantaneously  to  changes 
in  pressures  in  the  alveolar  chambers  such  that  no  net  change  in 
air  mass  occurs  within  the  airways.  This  assumption  is  justified 
because  the  pressure  front  due  to  a  change  in  pressure  at  any  point 
in  the  system  travels  to  other  points  at  the  speed  of  sound  although 
mass  transport  occurs  at  the  speed  of  flow  caused  by  the 
instantaneous  pressure  gradient.  Since  no  point  in  the  system  is 
much  over  1.5  feet  away  from  any  other  point  and  the  speed  of  sound 
is  1090  ft/sec  it  is  reasonable  to  assume  that  pressure  will  have 
redistributed  within  0.0015  sec.  If  the  time  step,  is  larger 
than  this  value  the  steady  state  assumption  will  be  not  only  adequate 
but  necessary.  The  modified  Newton-Raphson  iteration  described 
ear''ier  was  emoloyed  to  solve  the  steady  state  equations  for  the 
pressi.  es  within  the  airways  and  proved  to  be  very  satisfactory. 

Two  other  stability  criteria  must  be  imposed  on  the  integration 
P'ocedure,  namely; 

(10)  At-i^:, 

and 

(11)  At  ■'  RC  . 

Equation  OO)  simply  states  that  a  chamber  must  not  empty  before 
incoming  gas  has  had  a  chance  to  mix  with  resident  gas  and  Equation  (ll) 
is  just  the  resistance-compliance  criterion  imposed  uoon  any 
resi stor-rapaci tor  system.  These  conditions  allow  operation  of  the 
model  within  a  range  of  At  greater  than  0.002  sec.  for  forced  breathina 
and  less  than  0.025  sec.  for  passive  breathing. 


The  results  of  passive  breathing  simulation,  once  all 
parameter  values  have  been  set,  are  given  in  Figure  (7).  This 
figure  shows  simultaneous  values  for  average  alveolar  pressure. 
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as  '"Roresented  in  Figure  (7)  where  (P  ,  -P  ,  )  represents  the  ranae 

over  which  the  driving  pressure  will  operate,  t  is  the  elapsed  time 

during  the  breathing  phase  (inspiration  or  expiration),  t^^^  is 

the  phase  period,  and  y  is  a  constant  related  to  the  strength  of 

contraction  which  determines  the  rate  at  which  contraction  takes 

place.  and  P^.jp»  being  the  end  points  of  inspiration  and 

expiration,  may  be  initially  estimated  by  knowina  the  desired 

volume  chanoe  and  using  a  procedure  similar  to  that  used  to  calculate 

FRC,  assuminq  that  equilibrium  will  be  approached  and  all  gas  chamber 

pressures  will  be  near-  atmospheric.  If  equilibrium  will  not  be 

approached  during  the  run  P  and  P  .  will  have  to  be  stipulated 

max  min 

and  experimentally  determined. 


VOLUME  (L.)  FLOW  (  L  /  sec.  )  PRESSURE 


12  3  4  5 

TIME  (  s«Ci  ) 


FIGURE  7 


RESISTANCE  ( torry  L /sec.) 


30 


It  was  hoped  that  by  adjusting  y,  and  it  would 

be  possible  to  simulate  any  breathing  maneuver.  When  forced 
breathing  simulation  was  attempted,  however.  Equation  (12)  failed, 
yielding  a  peak  expiratory  flow  rate  far  too  low,  or  reaching  total 
lung  capacity  (TLC)  far  too  soon.  This  problem  will  be  discussed 
more  fully  in  the  next  two  chapters. 
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III.  RESISTANCE 

Classically,  airway  resistance  has  been  expressed  as 

(13)  AP  =  RV 

where  AP  is  the  pressure  difference  between  the  lungs  and  the 

mouth,  R  is  resistance,  and  V  is  volume  flow  rate.  This,  of 

course,  is  the  equation  for  laminar  flow  through  a  rigid  tube.  As 

( 39  40  41  ) 

early  as  1915,  however,  Rohrer  ’  ’  ,  who  was  the  first  to 

apply  fluid  dynamics  to  the  study  of  lung  mechanics,  noted  that 
pressure-flow  data  was  not  linear,  but  could  better  be  represented 
by  the  quadratic 

(14)  AP  =  KiV  +  K2V^ 

where  and  K2  ane  constants  characteristic  of  the  data  set. 

This  relation,  though  empirical  in  nature,  may  be  derived  from  fluid 

(42 ) 

mechanics  if  several  simplifying  assumptions  are  made  .  An 
analysis  of  this  derivation  brings  out  several  important  characteristics 
of  the  constants  K-j  and  K^;  but  first,  a  more  formal  expression 
of  the  Rohrer  equation  should  be  given.  Since  air  flow  may  be 
positive  or  negative  and  AP  must  always  reflect  that  sign. 

Equation  14  should  be  expressed 

(15)  AP  =  K^V  =  K2!ViV'  . 

Airway  resistance  now  becomes 

(16)  R  =  +  K2IVI  . 

(42 ) 

From  the  analysis  performed  by  Collins  ’  resistance  due 

to  laminar  flow  and  is  proportional  to  gas  viscosity  and  independent 
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of  density,  and  K2jV|  is  resistance  due  to  turbulent  flow  and 
K2  is  proportional  to  gas  density  and  independent  of  viscosity. 

Though  the  Rohrer  equation  is  often  cited,  it  is  seldom  used. 
Investigators,  in  general,  defer  to  the  laminar  flow  equation 
(Equaion  13)  and  assume  R  to  be  constant  for  the  sake  of  simolicity. 
Clinical  convention  would  have  resistance  recorded  when  flow  is 
near  0.5  liter/sec  though  this  is  not  rigidly  adhered  to  and  most 
often  flow  rate  is  not  mentioned  when  resistance  data  is  presented. 

Fry  and  Hyatt^^^  have  developed  a  model  to  study  the  mechanics 
of  airway  resistance  which  involves  a  much  deeper  analysis  than 
Rohrer's.  Their  model,  though  consisting  of  only  a  single  lung 
chamber  and  airpassage,  is  an  excellent  example  of  how  simple  models 
may  be  used  to  examine  and  explain  complicated  biological  orocesses. 
Their  model  (see  Figure  8)  emphasizes  the  elastic  character  of  the 
airoassaqes.  They  divided  the  conductive  zone  into  two  regions; 
those  interior  (region  I  in  Figure  8)  and  those  exterior  (region  E) 
to  the  lungs  proper.  During  expiration  the  airways  exterior  to  the 
lung  tissue  will  be  subject  to  high  pleural  pressures  on  the  outside 
and  low  pressures,  due  to  the  pressure  droo  down  the  airpassage. 
on  the  inside.  This  inward  directed  transmural  pressure  gradient 
will  cause  constriction  of  the  airpassage.  Since,  by  Poiseuille's 
Law,  resistance  is  inversely  proportional  to  the  fourth  power  of 
the  radius,  a  small  reduction  in  cross-section  of  the  tube  may 
increase  resistance  dramatically.  During  inspiration  the  opposite 
process  may  be  noticed. 

In  region  I  the  transmural  pressures  may  not  be  as  great  as 
those  in  reaion  E  as  region  I  is  subject  to  lung  pressure  rather 
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than  pleural  pressure;  but,  the  airpassages  in  this  region  are, 
in  general,  more  compliant.  Things  are  further  complicated  as 
these  airways  are  supported  oy  surrounding  iung  tissue  such  that 
their  diameters  will  be  affected,  not  only  by  their  own  elastance, 
but  also  by  pulmonary  elastance,  hence,  lung  volume. 

Results  from  the  Fry-Hyatt  model  indicate  that  for  a  given 
lung  volume  a  driving  pressure  snould  exist  which  yields  a  maximum 
flow  rate  (Figure  9).  If  the  driving  pressure  is  further  increased, 
flow  rate  will  either  plateau  or  decrease  due  to  airway  collapse. 
These  results  were  born  out  nicely  when  they  niade  measurements  on 
human  subjects  (see  Figure  10),  This  indicates  that  it  is  not 
enough  to  specify  the  flow  rate  at  which  a  tesistance  is  given;  the 
lung  volume  must  likewise  be  specified. 

While  Fry  and  Hyatt's  model  points  out  some  important  and 
interesting  aspects  of  airv.'ay  resistance,  it  is  not  easy  to  apply 
the.n  to  Rohrer's  resistance  equation  and  they  made  no  mention  of 
it.  T'np  Rohrer  constants,  constants  K-|  K2  in  Equations  (15) 

and  ll6)  cannot  be  experimentally  separated  such  that  data  describing 
either  independent  of  the  other  may  be  obtained.  As  an  attempt 
at  incorporating  Fry  and  Hyatt's  results  into  the  Oohrp>^  tnrmalism, 
this  author  used  the  followina  equations; 


17: 


B.|  + 
laB3V 


and 


1 3 : 


= 


I  *  By 


-•■Pa 


SO-VOLUME  PRESSURE  FLOW  CURVES 


where  -By  are  constants  and  P‘  is  the  transmural  pressure. 
These  equations  allow  for  an  increase  in  resistance  due  to 
increased  transmural  pressure  or  decreased  lung  \/olume.  The 
constants  B-j  -  By  were  set  experimentally  by  runiinq  the  model 
and  comparing  the  predicted  flow-volume  curves  obtained  with  those 
observed  cl inical ly. 

For  moderate  flow  rates  this  set  of  equations  performed 
adequately;  for  the  high  flow  rates  observed  in  the  maximum 
expiratory  flow-volume  (MEFV)  test,  however,  it  soon  became 
apparent  that  no  amount  of  adjusting  of  the  constants  B-j  -  By 
would  yield  the  desired  curve.  This  problem  was  partially  remedied 
by  modifying  the  shape  of  the  driving  function  from  that  of  a 
simple  exponential  rise  to  an  asymptote  represented  by  the  dashed 
line  in  Figure  (11)  to  a  function  yielding  the  torturous  curve 
represented  by  the  solid  line.  The  shape  of  this  final  curve  was 
justified  as  follows: 

a)  a  rapid  isometric  rise  in  pressure  to  near  maximal, 

b)  a  drop  in  pressure  as  inertia  is  overcome,  and 

c)  a  slow  isotonic  rise  to  maximum  pressure. 

This  was  disappointing  as  it  indicated  that  inertial  effects  of  the 
muscle  mass  which  were  assumed  small  in  developing  the  equations 
of  force  balance  were  indeed  significant  and  would  have  to  be 
incoroorated  into  the  driving  function  for  every  breathing  maneuver 
examined. 

While  attempting  to  model  emphysema  several  problems  surfaced 
which  caused  serious  doubt  as  to  the  validity  of  the  resistance  and 
compliance  equations.  To  model  a  mild  case  of  pulmonary  emphysema 
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in  the  absence  of  interstitial  fibrosis  it  should  only  be 
necessary  to  lower  the  compliance  of  the  affected  area.  When  this 
was  attempted,  however,  the  resultinq  FEFV  curves  did  not  display 
shapes  characteristic  of  emphysema  and  no  amount  of  adjusting  of 
the  compliance  curves  or  the  forcinn  function  seemed  to  help.  This 
forced  a  reevaluation  of  both  the  compliance  curves,  as  will  be 
discussed  in  the  next  chapter,  and  the  Rohrer  resistance  equations. 

At  this  point,  it  was  decided  that  rather  than  trying  to  pick 
new  equations  "out  of  hat"  to  replace  Equations  (17)  and  (18),  and 
adjust  them  to  match  flow-volume-pressure  data  with  no  guarantee 
of  success,  it  would  be  better  to  work  backwards  from  the  observed 
data  and  develope  an  equation  which  adequately  described  resistance. 
This  unfortunately  meant  the  Rohrer  resistance  equation  had  to  be 
abandoned  in  favor  of  a  curve  fit.  ihe  only  real  loss  entailed 
by  this  change  is  in  flexibility,  as  the  effects  of  altering  gas 
density  and  viscosity  cannot  be  studied  as  readily.  This  is  only 
a  minor  concession,  however. 

Fry's  data  (Figure  10)  is  given  in  three  dimensional  form  with 
coordinates  of  flow,  volume,  and  transmural  pressure.  This  can  be 
converted  to  resistance,  volume,  pressure  data  using  Equation  (13) 
where 
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P  =  P  ,  -  P  , 
pi  el 


St 


where  P^^  is  the  transmural  pressure  and  P 


el 


is  the  static 


St 


elastic  pressu’^e,  i.e.,  the  pressure  at  the  point  of  zero  flow  for 


a  given  volume.  When  this  transformation  was  performed  and  the 
data  plotted  lIic  family  of  curves  shown  in  Fiauren2)  was  obtained. 


AIRWAY  RESISTANCE  V.S.  ELASTIC  PRESSURE  FOR  DIFFERENT 

VOLUMES  ABOVE  R.V. 


O 
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■’’he  left  and  right  side  of  each  curve  may  be  described  as  the 
left  and  right  sides  of  a  pair  of  hyperbolas  that  share  a  common 
major  axis  on  the  resistance  axis,  a  common  center  and  focal  point, 
but  have  different  asymptotes.  The  two  sides,  representing 
inspiration  and  expiration,  of  each  curve  were  fit  to  the  hyper¬ 
bolic  function 

(20)  R=  E  f  M[-C-  + 


where  B  is  the  resistance  at  zero  flow,  C  is  half  the  minor 
axis,  and  M  is  the  slope  of  the  asymptote,  each  of  which  may  be 
apui oximated  graphically. 

Once  a  set  of  values  for  the  parameters  B,  C,  and  M  describing 
each  >.urve  were  obtained  the  volume  dependence  of  each  could  be 
modeled  by  a  curve  fit.  The  following  equations  were  chosen; 


(21) 

(22) 

(23) 

(24) 


B  =  6.472  -  1.665  V+  .2088  V' 
C  =  30.0 


M.n3p  =  0.0974/V 


exp 


,1 .48  ^  '  „2.5^ 


An  electrical  analog  (Figure  13)  v/as  employed  in  the 
distribution  of  these  parameters  to  the  individual  air  passages 
such  that  when  summed,  the  resistance  predicted  by  Equations  (20-24) 
was  obtained.  This  was  accomplished  by  making  a  series  of  simplifying 
assumptions.  As  a  first  aporoximation  the  lungs  were  assumed  to  be 
symmetric;  thus,  Figure  (13)  simpl  ifies  to  the  analogue  depicted  in 


3P//?0/^£re/? 


P//P  PPSSAGt'S 

- i 


ELECTRICAL  ANALOG 

figure  13 
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Figure  n4-A)where  the  subscripts  R  and  L  indicate  the  right 
and  left  sides  of  the  lung  and  is  equal  to  Rj^.  The  total 
resistance  of  the  system  is 


(25) 


R.  Rp 


^2  ^  ^ "  ^2  ?  *^1 


Substituting  Equation  (20)  into  Equation  (25)  wi  th  appropriate  subscripts 
yiel ds 


(26)  B+M[(C^+P^)^^^  -  C]  = 

B2+M2[(C2  +  -  C2]  ':B^+M^[(cf+P^)^^^  -  C^]; 


Since  P  on  the  left  of  the  equation  is  the  pressure  drop  throughout 
the  system, 

(27)  P  =  P2  +  Pl  . 

The  constant  C  has  units  of  pressure  and  it  was  assumed  that  when 
P  equals  C,  it  would  follow  that  P^^  would  equal  and  P^ 

would  equal  C2  such  that 

(28)  C  =  C2  +  Cl  . 

Further  examining  Equati on  (26),  we  find  that  when  P  equals  zero, 

(29)  B  =  B2+^Bl 
and  if  P  is  equal  to  C, 

(30)  M  =  ~  (MC  -  M2C2)  . 


Employment  of  Equations  (28-30)  assures  that  whatever  values  are 
assigned  to  the  distributed  constants  in  Equations  (21 -24),  for  the 
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resistance  R2,  the  total  resultant  resistance  will  be  approximately 
that  predicted  by  Equation  (20 J . 

The  equations  describing  the  distribution  of  resistance 
parameters  to  the  lower  regions  were  developed  in  a  similar  manner 
using  the  resistance  analogs  depicted  in  Figures  (14-B)  and  (14-C) 
and  will  not  herein  be  detailed.  It  is,  however,  important  to  list 
the  assumptions  made. 

1)  As  R2  represents  resistance  to  flow  through  air 
passages  exterior  to  the  lung,  it  should  not  be  a 
function  of  lung  volume. 

2)  The  dependence  of  resistance  on  volume  increases  as 
one  descends  into  the  lung  due  to  the  disappearance 
of  the  supporting  rings  of  hyaline  cartilage. 

3)  In  the  supine  position,  when  the  gravitational  pressure 
gradient  disappears,  the  total  pressure  drop  to  each 
alveolar  chamber  should  be  the  same. 

The  (Ciuioi  of  tiiic  cu:  «£  fit  are  illustrated  in  Figure  (12). 
where  the  digits  indicate  original  data^^^  and  the  solid  lines 
are  the  curve  fits.  The  effects  of  the  changes  in  the  handling  of 
airway  resistance  on  the  dynamics  of  t.f?  sy^tc-^..  p“^nrrially  the 
forcing  function,  shall  be  discussed  at  the  close  of  the  following 
chapter  on  compliance. 
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IV.  COMPLIANCE 


For  clinical  purposes,  compliance  is  defined  in  terms  of  a  sinqle 
chambered  lung  model  in  which  a  spring  represents  elastance^^'  For 
such  a  model.  Hook's  law  states  that  elastic  pressure  will  be  propor¬ 
tional  to  the  volume  displacement  from  the  equilibrium  volume,  V^, 
which  is  the  volume  the  lung  would  assume  if  elastance  were  the  only 
force  acting  on  it;  thus. 


(31) 


el, 


The  equation  describing  the  balance  of  forces  across  the  alveolar 
wall  in  such  a  model  is  just  Equation  (6)  here  repeated  as 


(32) 


al  V  '  pi  ' el 


where  the  hydrostatic  pressure  is  included  in  Pp.| .  The  equation 
desc'"'bing  flow  is  Equation  (13); 


(  *3  ^ 


where  is  pressure  at  the  mouth,  is  total  air  resiscance, 

V  is  volume  flow  rate,  and  inertance  is  assumed  negligible.  Combining 
these  three  eouations  yields 


(34) 


P 

mo 


+  R^V 


1 


(V-V^)  . 


Due  to  problems  in  determining  Ry,  compliance  is  always  determined 
at  instants  of  no  flow  where  Equation  (34)  simplifies  to 


P 

mo 


=  f-  (V-Vj)  . 

Cl  t 


(35) 
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Since  cannot  be  measured  indirectly,  zero  flow  must  be 

estab''  -.led  for  at  least  two  different  lung  volumes,  such  that 
may  be  determined  through  Equation  (3b)  as 


(36) 


where  is  assumed  constant  over  the  volume  range,  ,'.V.  This  equation 

is  just  the  difference  form  of  the  common  definition  of  compliance 
which  is 

(37)  E-i  ^  , 

where  E  is  elastance. 

The  clinical  procedure  for  measuring  using  Equation  (36) 
is  to  have  the  subject  inspire  a  small  volume  of  air  and  hold  it  until 
equilibrium  is  reached  throughout  the  lungs.  The  pleural  pressure  is 
measured  before  and  after  the  process  with  an  esophageal  balloon  (a 
slender,  very  compliant  balloon  inserted  through  the  nasal  passages, 
past  the  epiglotis,  and  positioned  midway  down  the  esophagus  which 
may  be  inflated  to  a  known  pressure  and  used  to  monitor  changes  in 
pressure  in  the  thoracic  cavity)  and  'J^  is  measured  with  a  spirometer. 
This,  after  a  number  of  such  measurements,  yields  the  static  pressure- 
volume  curve;  the  slope  of  which  is  "static"  compliance.  Figure  (15)^^^^ 
illustrates  such  a  curve. 

The  static  compliance  curve  indicates  clearly  defined  lung  volume 
maxima  or  total  lung  capacities  (TLCl  ot  between  5.5  and  7.0  liters, 
depending  on  the  subject's  height  and  age,  with  6.0  liters  being  the 
norm.  This  is  curiously  lower  than  results  obtained  frcni  excised 


o 


STATIC  PRESSURE  VOLUME  CURVE 
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lungs  by  many  investigators  as  reviewed  by  Hoppin  and  Hildebrandt^^^*^^ 

who  quote  figures  of  12  to  14  ml  gas  per  gram  lung  weight.  The  average 

f  45) 

lung  weighs  approximately  1.2  kg^  which  yields  an  expected  TLC  of 
14  to  17  liters,  more  than  twice  that  observed  clinically.  A  reasonabK 
explanation  for  this  difference  in  measured  TLC  has  not  been  found. 

During  the  early  scages  of  development  of  the  model  the  compliance 
curve  described  by  t.  e  exciseu  lung  experiment  was  selected  over  the 
invivo  clinical  measurement.  It  seemed  the  more  direct  procedure  and 
the  use  of  the  esophageal  Walloon  to  approximate  pleural  pressure  is 
subject  to  criticism  as  it  includes  the  compliance  of  the  esophageal 
wall.  Furthermore,  it  would  seem  that  lung  tissue,  being  more 
delicate  than  the  muscles  and  bone  surrounding  it,  should  not  furnish 
the  elastic  limits  of  the  system. 

The  equation  chosen  by  this  author  to  model  the  static  compliance 
as  depicted  in  Figure  (15)  is 


;38) 


P  . 
e  I 


tan  [ 


(V 


A2) 


2A 


+  A. 


where 


(39)  A  =  A^  -  A4 

when  the  volume  V  is  less  than  A2  and 

(40)  A=A3-A2 

for  V  greater  than  ^2',  A.j  is  the  slope  at  the  point  of  inflection, 
is  the  volume  at  the  point  of  inflection,  A^  is  the  lower 
asymptote,  A^  is  the  upper  asymptote,  and  A^ 


is  the  pressure  at 
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the  point  of  inflection.  The  compliance  constants  were  then 
distributed  to  the  various  pulmonary  elements  in  a  manner  analogous 
to  that  described  in  the  last  chapter  for  resistances.  All  lung 
■Jssue  chambers  are  assumed  to  be  identical  as  there  has  been  no 
evidence  indicating  different  characteristics  of  tissue  from  different 
regions  of  the  lung.  This  allows  that  in  the  supine  position,  where 
the  hydrostatic  pressures  are  equal,  all  lung  chambers  are  of  equal 
volume. 

As  stated  in  the  previous  chapter,  problems  arose  when  a  purely 
compliance  oriented  disease,  emphysema,  was  examined  which  necessitated 
the  reevaluation  of  the  compliance  philosophy.  This  was  accomplished 
by  developing  and  examining  a  simple  model  of  how  the  esophageal 
balloon  measures  pleural  pressure  (See  Figure  16).  Both  pleural  and 
esophageal  chambers  are  assumed  to  be  filled  with  incompressible  fluid 
and  the  balloon  is  filled  with  compressible  gas.  The  force  balance 
equations  for  the  two  pistons  are 


(41) 

and 

(42) 


P  ,  -  P  ,  +  P^ 
pi  el^  E 


B 


Combining  Equations  (41)  and  (42)  yields 


(43) 


Pr-P  1  -P.1  ■ 


pi  'el^ 


‘B 


Assuming  P  •,  is  to  be  negligibly  small,  a  good  assumption  considering 
B 

the  compliance  of  the  balloon  used,  and  reat ranging  Equf  .ion  (43)  yields 
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(44) 


B 


Since  the  balloon  is  inflated  slightly,  Pg-j  will  be  negative  in  the 
direction  of  the  arrow  in  the  figure.  This  implies  that  the  observed 
balloon  pressure  cannot  be  less  than  the  pleural  pressure  and, 
therefore,  the  esophageal  elastance  cannot  be  used  to  explain  the 
disparge  between  in-vivo  and  excised  lung  data.  Furthermore,  since 
the  upper  volume  limit  of  the  lung  cannot  be  less  than  that  indicated 
by  the  balloon  pressure  curve,  it  was  determined  that  esophageal 
compliance  must  be  relatively  large  such  that  P  ,  contributes 
little  to  Pg.  Pg  should,  therefore,  provide  a  good  approximation 
of  the  pleural  pressure. 

This  discovery  dictated  the  abandonment  of  excised  lung 
compliance  curves  in  favor  of  in-vivo  curves.  Figure  (17)  gives  the 
new  compliance  curve  generated  by  the  model  which  should  be  compared 
to  the  compliance  curve  for  a  human  subject  (See  Figure  15). 

The  initiation  of  the  modifications,  herein  descirbed  in  the 
compliance  equations  and  the  changes  in  the  resistance  equations 
discussed  in  the  previous  chapter  caused  one  unforseen  problem;  the 
modified  Newton-Raphson  method  of  solving  the  steady  state  pressure 
equations  in  the  airways  failed  to  converge  under  certain  conditions 
and  had  to  be  abandoned.  A  numerical  method  was  developed  whereby 
the  pressure  dependence  of  the  airway  resistance  is  neglected  during 
a  single  iteration.  This  allows  the  new  pressure  to  be  calculated 
from  the  mass-balance  equations  linearly.  Resistance  is  then 
calculated  in  terms  of  the  new  pressures  and  the  process  is  repeated. 
This  iteration  scheme  solved  the  problem  and  proved  more  stable  than 
the  modified  Newton-Raphson  method. 


PREDICTED  STATIC  PRESSURE  VOLUME  CURVE 
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The  solution  of  this  series  of  problems  proved  to  be  a  major 
advance  in  the  continued  development  of  the  model.  Not  only  did  it 
allow  the  examination  of  lung  mechanics  in  the  most  nonlinear  operational 
regions,  but  it  cleared  up  several  problem  parameter  areas  that  had 
nagged  at  the  credibility  of  the  model  since  the  first  runs.  As  has 
been  mentioned,  the  predicted  static  compliance  curve  now  provides 
a  much  better  match  to  clinical  data.  A  second  improvement  is  in  the 
formulation  of  the  driving  function  representing  contraction  of  the 
respiratory  muscles. 

As  was  discussed  in  the  previous  chapter,  the  examination  of 
forced  expiratory  breathing  maneuvers  has  been  hampered  by  the 
sensitivity  of  the  FEFV  curve  to  the  form  of  the  driving  functions. 

With  the  incorporation  of  the  new  resistance  and  compliance  changes 
this  problem  disappeared.  What  had  appeared  to  be  a  necessity  to 
incorporate  inertial  effects  of  the  muscle  masses  into  the  driving 
functions  proved  to  have  been  problems  with  the  resistance  and 
compliance  equations  that  could  be  partially  compensated  for  by  altering 
the  form  of  the  driving  function. 

Good  correlation  between  model  predictions  and  clinical  results 
for  a  number  of  forced  breathing  maneuvers  has  now  been  obtained,  this 
using  only  the  simple  exponential  form  for  the  driving  function  as 
used  for  passive  breathing  maneuvers.  Equation  (12)  with  appropriate  values 
for  the  constants.  Figure  (18)  illustrates  the  results  of  several 
simulations  of  forced  expiratory  breathing  maneuvers  using  different 
values  for  the  constant  y.  This  is  precisely  like  the  family  of 
curves  reported  first  by  Fry  and  Hyatt^^^^  for  a  series  of  graded 
effort  FEFV  maneuvers  (See  inset  on  Figure  18). 


VARIABLE  EFFORT  FLOW  VOLUME  CURVES 


FIGURE  18 


This  study  determines  values  for  two  important  lung  parameters 

which  need  setting  prior  to  the  WBA  studies.  These  are  the  maximum 

that  may  be  generated  by  muscular  contraction  and  the  maximum 

rate  oT  muscular  contraction,  i.e.,  P  ana  the  maximum  value  for 

max 

y.  Pj^^^  is  now  set  at  118  mm  Hg  and  the  maximum  value  for  y  is 
set  at  75. 


V.  Whole  Body  Acceleration 


As  stated  in  Chapter  I,  the  effects  of  WBA  on  the  human 

pulmonary  system  are  difficult  to  ascertain  using  traditional 

techniques  due  to  difficulties  in  collected  data  on  human  subjects 

and  due  to  the  inadequacy  of  animal  models.  It  is  proposed  that 

the  mathematical  simulator  of  the  human  pulmonary  system  already 

described,  if  supported  by  a  limited  amount  of  data  to  verify  the 

model's  preductions,  may  be  used  to  remedy  deficiencies  in  existing 

data  and  explain  observed,  but  poorly  understood,  phenomena  in 

terms  of  overall  lung  dynamics  as  well  as  regional  parameters 

heretofore  inaccessible.  In  this  chapter  we  will  attempt  to  use 

the  model  to  establish  a  basic  understanding  of  pulmonary  mechanics 

and  to  evaluate  a  technique  employed  to  reduce  the  adverse 

effects  of  WBA  and  improve  G-tolerance,  that  of  altering  the 

angle  of  attack  of  the  G-forces. 

It  has  long  been  recognized  that  at  some  seat  tilts  pilots 

are  able  to  withstand  greater  WBA  than  at  others.  U.S.  Air 

Force  studies  have  considered  angles  of  23",  28°  and  40°  as  the 

seatback  angles  for  the  F-15,  YF-16,  and  YF-17  fighter  planes. 

(47 ) 

Results  reported  by  Cell  indicate  that  seat  tilts  of  less  than 

45°  do  not  yield  significant  improvement  in  WBA  tolerance  and 
(48)  (49) 

both  Crossley'  '  and  Burns'  ’  report  that  for  seatback  angles 
of  65°  to  75°  tolerance  improves  sharply.  These  reports  are  for 
the  most  part  based  on  subjective  observations  of  pilots  and, 
therefore,  provide  only  very  loosely  defined  ranges  of  optimal 
seatback  angles. 
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In  the  attempt  to  gain  a  better  understanding  of  the  operation 
of  pulmonary  mechanics  under  WBA  and  how  changes  in  the  vector 
direction  of  the  applied  6-forces  alter  those  mechanics,  several 
lung  parameters  that  have  long  been  recognized  as  being  of 
importance  clinically  will  be  examined.  Lung  mechanics  may  be 
broken  down  into  two  broad  headings:  1)  Statics  and  2)  Dynamics. 
Under  the  topic  of  statics,  we  shall  examine  changes  in 

A)  Functional  Residual  Capacity  and 

B)  regional  distribution  of  resting  volume. 

Under  the  topic  of  dynamics,  we  shall  concentrate  on 

C)  Dynamic  Compliance  Loops  and 

D)  the  work  of  breathing. 

Study  A:  Functional  Residual  Capacity 

In  the  study  of  the  combined  effects  of  WBA  and  seatback  angle 
on  lung  parameters,  the  first  area  to  be  examined  will  be  that  of 
Functional  Residual  Capacity  (FRC).  This  is  the  equilibrium 
volume  of  the  lung  where  all  respiratory  muscles  are  in  their 
relaxed  state,  i.e.,  the  volume  of  the  lung  at  the  end  uf  relaxed 
expiration  from  which  active  inspiration  starts.  Consequently, 

FRC  defines  that  region  of  the  compliance  curve  over  which  the 
lung  operates  during  normal  breathing.  Since  normal  FRC  in  the 
erect  position  falls  near  the  most  compliant  portion  of  the  static 
volume-pressure  curve,  it  would  seem  that  maintainance  of  near 
normal  FRC  would  prove  most  conducive  to  WBA  tolerance. 


The  effect  of  change  in  body  tile  with  respect  to  eartn 
normal  gravity  was  once  an  area  of  intense  medical  interest, 
consequently,  much  data  is  available.  In  the  1920's  and  30':, 
it  was  well  known  that  changing  a  subject  from  an  upright  to  a 
prone  position  while  lying  in  a  bed  caused  a  reduction  in 
one  to  one  and  a  half  liters.  The  data  presented  as  plus 
symbols  in  Figure  09) are  taken  from  a  more  recent  source, 

Agustoni ,  but  are  exemplary  of  the  type  of  data  presentee  / 
early  researchers  and  illustrate  observed  changes  in  normal izec 
FRC  when  tilting  a  subject  from  an  upright  position.  O'  tilt, 
to  135'^  head  down  tilt.  Early  workers  capitalized  on  this  ana 
used  a  tilting  bed  to  respirate  subjects  with  paraiized  resoirator 
muscles.  Fortunately,  the  iron  lung  was  invented  and  this 
teennique  was  abandoned. 

The  dotted  line  in  Figure  09) i 1 lustrates  the  results  c* 
using  the  model  to  simulate  the  change  in  FRC  brougnt  about  ty 
body  tilt.  Clearly,  for  the  first  90''  to  100  moderately  gooa 
correspondence  with  the  literature  data  is  achieved.  Beyond  IOC 
tilt  the  model  predicts  that  FRC  will  increase  while  the  clinical 
data  continues  to  drop. 

Agosti  ni  ,  however,  demonstrated  that  FRC  measurement  is 
dependent  on  the  method  of  restraint  used  during  the  position 
changes.  He  points  out  that  the  "+"  symbols  are  characteristic 
of  data  collected  when  the  subject  is  supported  by  the  shoulders 
and  the  asterisks  are  characteri stic  of  data  collected  wnen 
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supporting  the  subject  by  the  ankles.  The  reasons  for  the  model's 
acting  as  if  it  were  supported  by  the  ankles  are  not  known  at 
this  time,  though,  it  must  be  a  characteristic  of  the  geometry  of 
the  components  chosen.  As  we  are  not,  at  this  time,  concerned 

with  bed  ridden  subjects  this  is  a  plus  rather  than  a  detriment. 

(51 1 

Glaister  reports  how  FRC  responds  to  changes  in  +G^ 
acceleration.  His  data,  including  error  bars,  are  reproduced  in 
Figure  (20).  Glaister  did  not,  however,  describe  the  details  of 
how  Lhe  subject  was  mounted  in  the  centrifuge.  It  may  be  assumed 
that  he  had  the  subject  seated  in  a  standard  pilot's  seat  with 
the  back  in  an  upright  position.  The  back  of  a  pilot's  seat,  as 
with  any  other  seat,  is  not  exactly  vertical,  but  is  set  at  some 
small  angle.  The  results  of  three  attempts  to  simulate  Glaister's 
results  using  different  seatback  angles  (0°  true  vertical,  5"', 
and  10")  are  presented  as  solid  lines  in  Figure  (20).  If  these 
assumptions  are  valid,  the  simulated  data  match  Glaister's  results 
nicely. 

If  the  model  has  the  capability  of  predicting  changes  in 
FRC  with  respect  to  changes  purely  in  body  angle  and  with  respect 
to  changes  purely  in  WBA,  the  model  provides  a  valid  tool  for 
the  examination  of  the  combined  effects  of  seatback  angle  and  WBA 
on  FRC.  Figure  (21)  is  a  three  dimensional  representation  of  such 
a  study  where  seat  tilts  of  0"  to  90"  and  WBA  of  from  1  G 
earth  normal  to  10  G's  are  examined.  The  1  G  iso-gravity  curve 
is  the  same  as  that  presented  in  Figure  (19)  and  the  0"  to  10" 
iso-angle  curves  are  the  corresponding  curves  from  Figure  (20). 


M0RMALI2ED  FUC  <  FRC  /  FRC  Ct  If.  0* ) 
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If  we  accept  our  earlier  hypothesis  that  maintaining  near 
normal  FRC  while  undergoing  high  sustained  WBA  is  ideal,  this' 
study  indicates  that  any  seatback  angle  from  15"’  to  65°  is 
acceptable.  Figure  (21),  however,  provides  only  gross  lung  data 
and  before  a  "best  angle"  is  chosen,  the  behavior  of  individual 
chambers  should  be  examined. 

Study  B:  Regional  Distribution  of  Resting  Volume 

Figure  P2)is  a  breakdown  of  the  data  from  Figure  (21)  into 
regional  values.  Each  region  has  a  maximum  volume  of  2.28  liters 
and  a  minimum  volume  of  0.22  liters.  It  may  be  noted  that  for 
1  G  and  0°  tilt,  the  upper  chambers  have  more  than  tw^'.e  the 
volume  of  the  lower  chambers.  This  is  a  reflection  of  the  increase 
in  hydrostatic  pressure  in  the  pleural  fluid  caused  by  the  1 
force  as  was  first  measured  by  Glazier'  ^  (see  page  15). 

Glazier  further  demonstrated  that  when  G^  was  increased,  this 
effect  was  magnified.  This  fact  corroborates  the  niodels  results 
when  the  0°  iso-angle  curves  for  each  region  are  compared. 

Examining  the  predictions  for  the  lower  chambers,  we  find 
that  increased  WBA  has  caused  a  reduction  in  regional  volume  to 
near  the  minimum  volume.  This  would  place  this  region  on  the 
lower  end  of  the  pressure-volume  curve  where  it  would  be  very 
noncompl iant.  Changes  in  seatback  angle  do  not  have  appreciable 
effect  in  this  region.  The  upper  chambers  behave  much  differently 
in  that  for  low  seat  tilts,  increase  in  WBA  causes  the  volume  to 
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LOCAL  FRC  IN  THREE  LUNG  REGIONS 

FIGURE  22 
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approach  its  maximal  value.  This  trend  is  not  appreciably  affected 
until  seatback  angle  reaches  a  very  narrow  range  of  between  60” 
and  70°  where  FRC  begins  to  plunge  toward  its  lower  limit.  At 
angles  of  80°  to  90°,  the  upper  chamber  is  also  near  its  minimum 
volume.  The  middle  chambers  offer  the  most  interesting  study  in 
that  there  are  two  possible  optimum  angles:  10°  and  60°. 

In  summation,  the  gross  lung  FRC  data  presented  in  Figure  (21) 
do  not  provide  the  detail  necessary  to  determine  an  optimum 
seatback  angle.  It  does,  however,  provide  results  which  may  be 
used  to  confirm  the  ability  of  the  model  to  simulate  static  lung 
mechanics.  An  examination  of  regional  FRC  provides  two  angles  as 
possible  candidates  for  consideration  as  optimum  seatback  angles, 
10°  and  60°  to  70°,  i .e. , 

i)  if,  due  to  other  considerations,  it  seems  advisable  to 
force  the  best  ventilation  to  occur  in  the  middle 
chambers,  this  may  be  accomplished  by  choosing  an 
angle  near  10°,  or 

ii)  if  it  be  deemed  advisable  to  share  the  best  ventilation 
between  all  of  the  chambers,  then  an  angle  of  between 
60°  and  75°  should  be  considered. 

A  choice  between  these  alternatives  cannot  be  made  on  the 
basis  of  a  single  study  but  will  be  reconsidered  in  the  next 
chapter  when  more  information  is  available. 
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Study  C:  Dynamic  Compliance 

In  Chapter  IV  the  concept  of  static  lung  compliance  is 
developed  and  the  measurement  technique  is  discussed.  As  was 
therein  mentioned,  compliance  is  determined  at  instants  of  zero 
flow  such  that  resistance  is  eliminated  from  the  measurement. 

A  second  technique  for  approximating  the  compliance  utilizes 
the  instant  of  zero  flow  at  the  mouth  that  occurs  naturally  at 
the  beginning  and  end  of  inspiration.  This  procedure,  which 
determines  the  "dynamic  compliance,"  calls  for  the  subject  to 
breath  small  tidal  volumes  around  FRC,  such  that  compliance 
should  change  little  throughout  the  breath.  The  procedure  may, 
however,  be  criticized  as  the  value  obtained  is  not  pure 
compliance.  The  fact  that  flow  has  ceased  at  the  mouth  is  no 
indication  that  flow  has  ceased  throughout  the  lungs.  Since  lung 
volume  and  pleural  pressure  are  out  of  phase  with  each  other  and 
since  dynamic  compliance  is  determined  during  a  dynamics  rather 
than  a  static  state,  flow  will  be  occurring  in  the  peripheral 
airways  (those  nearest  the  acini)  and  Equation  (34)  rather  than(35) 
should  be  employed  to  determine  Cj^. 

Though  this  criticism  indicates  that  the  value  determined 
for  is  not  a  compliance  in  the  true  sense,  it  leaves  us 

with  some  interesting  consequences.  Flow  is  occurring  and  dynamic 
compliance  will,  therefore,  be  a  function  of  the  resistance  and 
the  flow  rate.  Since  the  procedure  dictates  that  flow  in  the 
upper  airways  be  near  zero,  a  change  in 


must  be  a  function 
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of  peripheral  airway  resistance  only.  Peripheral  airway  resistance 
must  change  drastically  before  the  total  resistance  is  significantly 
affected,  but,  by  comparing  dynamic  compliance  to  values  accepted 
as  normal,  abnormal  peripheral  airway  resistance  may  be  detected 
in  diseases,  qualitatively,  far  earlier,  thus  giving  advance 
warning  of  the  onset  of  diseases  such  as  asthma  and  chronic 
bronchitis  which  seem  to  attack  primarily  small  airways.  The  use 
of  dynamic  compliance  in  determining  peripheral  airway  resistance 
is  further  refined  by  the  recognition  that  the  phase  shift  between 
and  Pp^  increases  with  the  frequency  of  breathing.  This 
means  that  the  flow  in  the  small  airways  will  increase  and 
dynamic  compliance  will  better  reflect  abnormal  resistances  in 
these  regions. 

Figure  (23-a)  i 1 lustrates  the  results  of  a  three  breath  run 
starting  at  FRC  in  which  the  resistance  was  normal  and  breathing 
frequency  was  30  BPM  (breaths  per  minute).  Figure  (23-b)  shows 
a  similar  run  in  which  the  linear  resistance  in  the  lower  order 
bronchi  feeding  the  lower  two  thirds  of  the  lung  was  increased 
tenfold.  The  drop  in  compliance  and  the  widening  of  the  P-V 
loop  expected  clinically  as  a  result  of  high  small  airway 
resistance  is  evident.  The  increase  in  total  resistance  associatea 
with  this  large  increase  in  lower  order  resistance  was  only  by  a 
factor  of  two,  however,  and  remained  within  the  normal  range. 

Such  an  abnormality  might  go  unnoticed  if  further  tests  were  not 
run  on  the  subject.  A  study  of  frequency  dependence  of 
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performed  on  the  same  subject,  provided  in  Figure  (24),  more  clearly 

indicates  the  problem.  The  shape  of  the  high  resistance  curve 

very  nicely  matches  the  experimental  results  observed  by  Ingram, 

et.al  i“)in  their  study  of  apparently  healthy  smokers  which  are 

included  as  an  inset  to  Figure  (24]l  The  shaded  area  in  the 

inset  is  the  range  in  which  compliance  of  the  normal  non-smokers 

fell.  (Note:  The  study  using  the  model  only  included  frequencies 

up  to  80  BPM  and  is  shifted  up  with  respect  to  the  inset  values 

due  to  a  different  definition  of  normalization,  i.e.,  C„  /Cp, 

Dyn  Dyn2Q 

rather  than  ) 

Due  to  the  difficulty  of  obtaining  dynamic  compliance  loops 
for  subjects  undergoing  WBA,  no  information  has  been  published 
indicating  how  WBA  affects  dynamic  compliance.  If,  however,  we 
accept  that  the  model  is  adequately  simulating  system  dynamics 
under  earth  normal  gravity  as  indicated  by  the  above  mentioned 
dynamic  compliance  studies,  the  passive  breathing  study  (Chapter  II), 
and  the  variable  effort  forced  maximal  expiration  study  (Chapter  IV), 
and  if  we  accept  that  system  statics  are  being  adequately  simulated, 
as  indicated  by  the  arguments  presented  in  Studies  A  and  B  of 
this  chapter,  the  model  should  provide  a  good  simulation  of  system 
dynamics  under  WBA. 

In  Figure  (25) are  given  the  results  of  a  series  of  simulated 
passive  breaths  while  in  different  body  positions  ranging  from  0° 
(upright)  to  90°  (supine)  and  for  different  levels  of  WBA  ranging 
from  1  G  to  7  G's.  Passive  breathing  was  defined  as  a  two 
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second  inspiration  of  about  0.5  liters  of  air  starting  at  FRC 
and  a  3  second  relaxed  expiration.  Each  loop  represents  a  pressure 
volume  plot  with  each  tic  mark  increasing  horizontally  representing 
2  mm  Hg  and  each  tic  mark  descending  vertically  representing 
0.1  liters  of  air  inspired.  (Inspiration  is  recorded  as  a  negative 
change  in  spirometer  volume.) 

It  is  interesting  to  note  that,  though  WBA  tends  to  decrease 
dynamic  compliance,  the  two  angles  most  affected  are  0°  and  90° 
with  90°  by  far  the  worse.  This  result  was  anticipated  from 
the  results  of  the  local  FRC  study  described  in  Study  B  of  this 
chapter.  The  FRC  of  each  chamber  is  changing  and  the  region  is 
constrained  to  operate  on  a  different  portion  of  its  static 
compliance  curve. 

Not  only  does  increased  WBA  cause  an  increase  in  dynamic 
compliance,  it  also  causes  changes  in  the  overall  shape  of  the 
pressure-volume  loop.  Shapes  such  as  those  suggested  for  higher 
WBA  and  seatback  angles  of  30°  to  60°  have  neven  been  observed 
before.  These  results  should  be  expected,  however,  as  the  lung 
is  no  longer  operating  on  the  most  linear  portion  of  the  static 
compliance  curve. 

Remembering  that  the  width  of  the  loop  is  an  indication  of 
overall  airway  resistance  we  may  evaluate  this  series  of  runs  in 
terms  of  both  compliance  and  resistance;  the  two  major  characteristics 
of  lung  tissue  mechanics.  Where  dynamic  compliance  is  least 
affected  by  WBA  for  tilts  of  15°  to  75°,  resistance  appears  least 
affected  for  tilts  of  0°  to  15°  and  60°  to  75°. 
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The  net  results  of  these  observations  match  those  ootained 
by  casual  visual  inspection  of  Figure  (25),  i.e.,  two  seatback 
angles  stand  out  as  being  optimal  for  minimizing  the  effects  of 
WBA  on  the  dynamic  compliance  angle,  the  width,  and  the  overall 
shape  of  the  pressure-volume  loops:  15°  and  75°. 

Study  D:  Work  of  Breathing 

Since  the  concern  of  pulmonary  mechanics  is  the  movement  of 
air  into  and  out  of  the  lungs,  the  mechanical  work  involved  may 
ta  expressed  as 

(45)  W  =  /  P  dV  . 

To  obtain  the  true  work  of  breathing,  P  should  be  the  driving 
pressure  generated  by  the  respiratory  muscles.  This  pressure 
cannot  be  measured,  however,  and  the  pressure  differe.ice  between 
pleural  pressure  and  mouth  pressure  as  measured  by  an  esophageal 
balloon  is  conventionally  used.  This  would  represent  the  area 
under  the  volume-pressure  curve  of  the  dynamic  compliance  loop 
and  would  only  yield  the  work  performed  in  overcoming  lung  tissue 
elastance  and  airway  resistance,  which  is  of  only  limited  value. 
The  pressures  involved  in  driving  the  model  are  known  and, 
consequently,  the  work  of  breathing  may  be  calculated  as 

(46)  «  =  “''tx  "  ^'’OA  “''da  ■ 

In  general,  when  biological  systems  are  being  examined,  it 
is  not  mechanical  work  but  expended  energy  or  effort,  of  which 
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mechanical  work  is  but  a  part,  which  is  of  interest.  Effort, 
unfortunately,  is  both  difficult  to  define  and  more  difficult 
to  measure  in  the  invivo  state.  If  the  entire  organism  is  being 
examined,  O2  consumption  and  heat  generation  may  be  used  as 
indicators.  When  only  certain  muscles  of  the  organism  are  of 
interest  and  especially  when  that  organism  is  man,  who  objects 
to  invasive  measurements,  expended  effort  is  virtually  impossible 
to  define  and  mechanical  work  is  often  used  as  an  indicator. 

Investigators  who  report  the  work  of  breathing  often  ignore 
these  facts  and  a  few  simple  examples  may  be  in  order  in 
clarifying  the  difference  between  the  two.  Were  the  equilibrium 
volume  of  either  the  diaphragm  or  the  chest  to  be  forced  either 
up  or  down  the  compliance  curve  to  near  the  asymptote  (see 
Figure  (15),  it  would  require  a  tremendous  effort  to  inspire  even 
a  small  volume  with  the  muscles  of  that  anatomical  element.  The 
work  performed  by  those  muscles  would,  however,  by  Equation  (46), 
be  small.  If,  on  the  other  hand,  the  system  were  to  be  driven 
entirely  by  one  muscle  complex  and  the  other  were  to  remain  passive 
and  exoend  no  effort,  the  pressures  in  the  system  would  cause 
movement  of  the  other  and  Equation  (46)  would  indicate  that  work  has 
been  performed  though  no  effort  had  been  expended. 

Figure  (26)  gives  the  results  of  the  calculated  work  of 
inspiration  for  the  simulated  breathing  studies  reported  in 
Figure  (2^.  Since,  for  relaxed  breathing,  expiration  is  accomplished 
by  passive  relaxation  of  the  respiratory  muscles  and  little  effort 


Work  of  Potsivo  Breathing 


figure  26 


is  actually  expended,  though  mechanical  work  near  that  of 
inspiration  is  performed,  only  the  work  of  inspiration  will  be 
examined.  To  make  comparison  to  data  easier,  both  muscle  groups 
were  constrained  to  exert  equal  driving  pressures.  Work  per 
liter  of  inspired  air  is  reported  as  not  all  breaths  were  of 
exactly  the  same  depth.  Ignoring  the  criticisms  of  using  work 
as  a  measure  of  effort  cited  earlier,  it  is  easy  to  interpret 
this  family  of  iso-G  curves.  Two  ranges  of  seatback  angles 
appear  best  in  minimizing  the  work  of  breathing;  best  is  the 
70°  to  75°  range,  and  surprisingly,  of  secondary  importance  is 
the  0°  to  15°  range. 

Two  points  are  of  special  interest  in  examining  Figure  (26) 
First,  note  how  independent  of  seatback  angle  the  work  of 
breathing  is  in  the  earth  normal  1  G  environment.  It  would 
seem  the  system  was  designed  to  operate  under  this  condition. 
Secondly,  it  is  interesting  to  note  how  poorly  the  system 
responds  to  WBA  in  the  supine  position.  This  position  has  been 
chosen  by  NASA  for  the  astronauts  primarily  on  the  basis  of 
circulatory  considerations. 


VI.  Conclusion 


In  the  preceding  chapter  we  have  described  four  studies  of 
the  effects  of  WBA  on  pulmonary  mechanics  using  the  mathematical 
simulator  developed  and  justified  in  the  preceding  chapters. 

Each  study,  to  the  extent  possible,  is  divided  into  two  parts; 

a)  verification  of  simulation  results  by  comparison  with 
existing  experimental  data  and 

b)  prediction  of  the  effects  of  WBA  on  lung  mechanics  in 
areas  hitherto  unexamined. 

These  will  now  be  discussed  in  more  general  terms  as  they  relate 
to  each  other. 

A.  Verification  and  Criticism  of  Results: 

While  the  first  of  these  studies,  that  of  the  effect  of  WBA 
on  gross  lung  FRC,  provides  the  least  useful  data  in  terms  of 
defining  optimum  seat  tilt,  it  is  perhaps  the  most  significant 
as  it  is  the  most  directly  verifiable  and  provides  the  best 
indication  that  the  model  is  adequately  simulating  the  system. 

The  second  study,  that  of  the  effects  of  WBA  on  local  equilibrium 
volumies,  provides  a  basis  for  evaluating  and  understanding  other 
studies.  While  it  is  supported  by  studies  involving  frozen  dogs, 
it  derives  its  main  support  from  the  fact  that  it  contributes 
the  parts  from  which  the  previous  study  involving  the  gross  TRC 
is  compiled. 
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As  we  are  primarily  interested  in  hw  well  the  lung  can 
function,  the  last  two  studies,  those  in  Iving  dynamics,  must  be 
given  the  greatest  weight.  Study  C  of  the  last  chapter  indicates 
that  good  correlation  between  model  results  and  experimental 
data  have  been  achieved  for  the  earth  normal  gravity  studies; 
however,  due  to  the  difficulties  involved,  no  data  are  available 
to  corroborate  the  high  WBA  studies.  Likewise,  no  data  are 
available  to  corroborate  the  work  of  breathing  studies  and 
Study  D  of  the  previous  chapter  discussed  many  of  the  pitfalls 
of  using  these  results  in  defining  difficulty  in  breathing. 

B.  Results 

Recognizing  the  above  criticisms  of  the  limitations  of  the 
studies  thus  far  carried  out,  it  is,  however,  necessary  to  make 
a  final  conclusion  based  upon  the  data  available  for  selecting 
an  optimum  seatback  angle  for  maintaining  good  respiratory  function 
under  WBA.  Summarizing: 

Study  B  indicates  a  choice  between  10°,  for  which  ventilation 
will  be  better  in  the  middle  chambers,  and  an  angle 
between  60°  and  70°  for  which  ventilation  would  be 
more  equally  distributed. 

Study  C  suggests  that  angles  of  15°  and  75°  both  provide  a 
near  normal  dynamic  compliance  loop. 

Study  D  predicts  that  an  angle  of  near  70°  is  best,  but 

that  an  angle  of  between  0°  and  15°  is  almost  as  good. 


with  respect  to  minimizing  the  work  of  breathing. 
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Clearly,  two  angle  ranges  have  emerged  as  candidates  as  the 
optimal  angle:  0°  to  15°  and  60°  to  75°.  No  clear  choice 
between  the  two  can  be  made  based  purely  on  the  pulmonary 
function  studies  performed  thus  far  and  other  factors  must  be 
considered. 

The  primary  cause  of  pilot  failure  is  blackout  which 

occurs  as  a  result  of  inadequate  delivery  of  O2  to  the  brain. 

This  may  be  the  result  of  inadequacy  of  one  or  more  of  the 

following  physiologic  phenomenon:  1)  pulmonary  ventilation, 

2)  blood-gas  communication,  and  3)  circulation.  The  first  of 

these  is  the  subject  of  this  dissertation  and  has  brought  us 

to  the  current  dilerma.  As  to  the  second,  it  is  known  that 

increased  WBA  affects  the  blood  and  gas  portions  of  the  pulmonary 

system  as  opposites,  i.e.,  the  increased  hydrostatic  gradient  in 

the  circulatory  system  causes  pooling  of  blood  in  the  lower 

regions  of  the  lung  while,  at  the  same  time,  it  causes  pooling  of 

gas  in  the  upper  regions,  as  discussed  in  Study  B.  This  same 

study  indicates  that  while  the  gas  is  pooling  in  the  upper 

regions  in  an  upright  posture,  the  best  ventilation  is  occurring 

in  the  middle  region.  This  result  is  interesting  when  compared 

(52) 

to  perfusion  data  obtained  by  Whinnery'  '  who  introduced  radio¬ 
active  microspheres  into  the  blood  stream  of  miniature  swine 
undergoing  high  sustained  +G^  acceleration.  The  microspheres 
lodged  in  the  pulmonary  capillary  bed  and  scintiscanning  was  used 
to  obtain  a  "picture"  of  pulmonary  perfusion.  His  results  indicate 
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that  the  best  perfusion  during  WBA  levels  of  2-6  G's  occurs 
about  half  to  two-thirds  of  the  way  down  the  lung.  This  is  the 
same  area  that  our  model  predicts  to  be  the  best  ventilated  under 
+G^  acceleration. 

Though  no  data  exist  in  confirmation,  it  would  seem  logical 
that  as  the  angle  of  attack  of  the  WBA  increases  toward  the 
supine  position  and  the  hydrostatic  pressure  difference  between 
top  and  bottom  of  the  lung  decreases,  perfusion  would  better 
distribute  blood  through  the  lung  just  as  ventilation  does.  We 
are  left  with  the  conclusion  that  an  examination  of  the  blood  gas 
communication  does  not  help  very  much  in  choosing  between  the  two 
seat  tilt  ranges  as  it  seems  that  the  location  of  the  best 
ventilation  and  perfusion  will  match  for  any  seat  tilt.  It  would 
appear,  on  the  surface,  that  the  60°  to  75°  angle  would  be  the 
best  as  it  allows  both  ventilation  and  perfusion  to  occur  over  a 
broader  region,  but  this  has  not  been  shown. 

The  last  of  the  physiological  phenomena  which  may  help  in 
choosing  the  best  angle  is  that  of  circulation.  This  is  the  area 
that  has  been  most  intensely  studied  and  all  data  points  to  the 
same  logical  conclusion:  in  the  upright  position,  increased  WBA 
causes  blood  to  pool  in  the  lower  extremities  and  inadequate  blood 
flow  to  the  brain  results.  Further,  an  increase  in  seatback  angle 
toward  the  prone  position  relieves  the  problem.  Taking  this  into 
consideration  we  must  conclude  that  the  60°  to  75°  seat  tilt 
should  be  optimum  in  relieving  the  effects  of  high  sustained  WBA. 
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The  studies  herein  discussed  form  a  basis  for  what  is  hoped 
win  be  a  long  and  detailed  examination  of  the  high  WBA 
environment.  In  the  process  of  determining  what  may  be  the 
optimum  seatback  angle  for  relieving  the  effects  of  WBA,  we  have 
developed  a  data  L’=i=  desc!  itirij  h«.»y  the  lung  operates  du;  .ng  pas-ive 
breathing.  This  may  be  used  for  comparison  for  future  studies. 

The  model  is  fully  capable  of  simulating  the  effects  of  various 
anti-G  suits  and  of  different  breathing  maneuvers  such  as  positive 
pressure  breathing  or  the  valsalva  maneuver.  It  is  also  capable 
of  simulating  the  effects  of  WBA  on  different  body  types  such  as 
tall  or  short  subjects,  slender  and  husky  builds,  and  obesity. 

A  long  range  goal  may  be  to  define  an  optimum  subject  type 
for  withstanding  WBA. 

The  model  is  not  limited  to  the  study  of  the  high  G 
environment.  As  the  study  of  the  effects  of  high  airway  resistance 
indicates,  this  model  should  be  very  useful  in  studying  various 
lung  diseases  and  developing  clinical  tests  for  their  early 
detection  and  diagnosis. 
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APPENDIX  A:  Derivation  of  Equations 


PART  I:  Alveolar  Equations  -  Table  I 

Restating  the  force  balance  equations  (  Equations  3,  4,  and  5  in 
the  text  )  we  have 

(A-1)  ^  ^  ^  ^ 

and 


(A-3)  p.  -  Pp/  f  y  P^/^  . 

Subtracting  Equation  A-3  from  A-1  to  eliminate  P^^^  yields 

(A-4)  P  ^  ^  P^,.  =  ^  Pg  • 

The  resulting  equation  has  been  rearranged  such  that  all  user  defined 
or  constant  functions  have  been  grouped  on  the  right.  These  we  rename 


(A- 5) 


''sj  ~  Ptk 


p  -b  P  -  P 


Differentiating  Equation  A-4  with  respect  to  time,  we  have 


(A-6) 


P 

J 


a^r. 


<aPe/ 


'J 


0 


(  Unless  otherwise  specified,  differentiation  will  always  be  assumed 
to  be  with  respect  to  time  which  will  be  Indicated  by  the  dot  notation.) 
Differentiating  the  equation  of  state  for  each  alveolar  chamber  j 
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yields 


U-,,  ^  4  . 

V 

Solving  f?**  Vj  from  EousClon*  A-6  erd  A-7  yield* 


0 

p-  -- 

7~ 

V- 

1/ 

a  y^rx 

• 

J 

will  now  be  eliminated  from  Equations  A-1  and  A-2  such  that 

(A-9)  ■  ^Tx  ~  ^a/t  ' 

The  user  defined  functions  are  again  on  the  right  and  will  be  renamed 


(A-IO)  --  ^  ^  ^  . 


Differentiating  Equation  A-9  yields 


(A-11)  di/r,  '  ay  DA  ~ 


SZ 


Differentiating  Equation  2  in  the  text,  the  volume  equation,  gives 
(A-12)  --  /T  ;>■  -  . 

m  • 

Eliminating  between  Equations  A-11  and  A-12  and  solving  for 


r  i-  Ps2, 

^  j 


(A-13)  i/, 


Tx 


^  DA  ^ 

The  equation  for  mass  balance  of  an  individual  gas  in  the  alveoli 


is  a  modified  form  of  Equations  11  and  12  in  the  text  such  that 
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(A- 14) 


there  the  subscript  m  Identifies  the  chamber,  i  or  j  ,  from  which 
mass  is  flowing  and  is  dependent  on  the  sign  of  .  Differentiating 
the  equation  of  state  for  gas  1  we  have 


(A-15)  Ci-  --  if'i,  ■  Cj, 


As  was  stated  in  the  text,  water  vapor  la  assumed  to  be  at  the 


saturation  point  such  that  if  the  temperature  is  constant,  C^j  is 

(46) 

a  function  of  P  only,  i.e.. 


(A-16) 


V 


^7-  ~  1 


'J 


^ere  is  the  molar  voltsne  of  water,  is  a  reference  pressure 

and  Pg  is  the  vapor  pressure  of  water  at  Pq  .  Since 


(A-17)  --  Cu.^-  l/j 

it  follows  that 

(A-18)  --  l/^-  P^-  ^ 

# 

As  n^j  is  not  just  a  function  of  mass  transfer  between  chambers, 
a  source  equation  representing  evaporation  and  condensation  oust 
be  supplied: 

4 

(A-19)  //  -  P- 


with  the  fact  In  mind  that  moles  of  water  must  be  treated  differ¬ 
ently  than  the  moles  of  other  gasses,  we  restate  the  equation  of  state 
as 

/P  /  / 

(A-20)  ^  pr  (  ^  ^7- .  y 

where  we  define  n^j  as 
(A-21)  IT  ■ 

Differentiating  Equation  A-20  yields 

(A-22)  /i.  -  ^ -  y  .  cjj/ 

where 

- 

(A-23)  ~  // 

V  ^ 

Equations  A-8,  A-13,  A-14,  A-15,  A-18,  A- 19,  and  A-22  esc=7;pass 
the  equations  listed  In  Table  I  of  the  text.  This  set  >f  equations 
Is  severely  algebraically  coupled.  By  algebraically  coupled  It  Is 
meant  that  If  each  unknown  (  the  set  of  symbols  on  the  left  side  of 
the  equations  In  Table  1  )  Is  treated  as  an  algebraic  variable  rather 
than  a  differential  function,  a  set  of  n  equations  In  n  unknowns 
Is  obtained.  This  system  of  equations  must  be  algebraically  manlp- 
tilated  until  It  can  be  solved  sequentially.  This  was  accomplished 
as  follows. 

Eliminating  P.  from  equations  A-18  and  A-22  yields 
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(A-24)  -  [  X^-  ^  ^  C 


where  Che  following  IdenClCiea  are  used 


-xc^)J 


7  / 


-  _ 

(A-25)  X,  .  XX  r  — . 


(A-26)  C,  --  X 


Further,  defining 


(A-27) 


Cu^-  -  X^ 


Equaclon  A-24  becomes 


(A-28) 


Allowing  Che  following  Identities, 


<A-2,)  xx,  -  I  , 


(A-30)  /9^ 


y  ~x 

i-  i3ji 

^  3  Xj 


Equation  A-8  becomes 


(A-32)  ^  ^  (f  >^4.^  ^  J  ^  ^ 
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Combining  Equations  A-28  and  A'~32  to  eliminate  n^j 


(A-33)  ^  ^ 


~trx- 

/  -  py 

/-Xj 


Applying  the  Identities 


(A-34)  ^ 


and 


(A-35)  CJ^- 


/- 


Xi£l 


IX j  r>r.- 
/  -  ‘^1 


/-Xy 

to  eqtiatlon  A>33,  we  have  the  much  simpler  equation 


(A-36)  7‘  OJj 


Finally,  substituting  the  Identities 

^ P/! 


(A-37)  j; 


and 


y. 


(A-38)  ^ 

Into  equations  A-13  gives  us 


(A-39)  K75, 


=  t  H  ^  * 


r-/ 


Combining  Equations  A-36  and  A-39  to  eliminate  Vj 


yields 


yields 
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,  ^  6Jr  ^  ^ 

(A-AO)  --  - - -  ^ 

/  -  t  O' ^ 

where  we  make  use  of  the  Identities 


(A-41) 


a  nd 


(A-42)  --  ZT  //  . 

The  ?y«»-«jn  of  equations  listed  In  Table  I  may  now  be  solved 
using  the  set  of  sequential  equations  listed  In  Table  IV. 


ftLl^EOL/lR  EGfUAT  IONS 


/. 

/9ji  -  r 

V 

Ej  - 

- 

it  OJj-  f-  ^ 

7  ~  <f-  y^- 

J 

■i. 

V. 

- 

/j-  = 

^  1  "T  '  ’  ‘^’  ^  ^  ^ 

4-  ■• 

{ - )  f 

7. 

c^  -  -■ 

Y 

X 

- 7  z^- 

f. 

‘H  ' 

table  iz 
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PART  II:  Airway  Equatloaa  -  Table  II 

Since  volumes  in  the  airway  chambers  are  assumed  constant,  the 
derivative  of  the  equation  of  state  for  gas  Is 


(A-43)  P 


iy. 


J 


)  • 


The  mass  balance  equation  Is  just  Equation  12  In  the  text.  The  deriv¬ 
ative  of  the  concentration  equation  for  gas  1  and  water  vapor  w 
will  be 


(A-4A)  £ 


(A-45) 


P- 


The  molar  values  of  water  In  each  chamber  will  then  be  defined  by  the 
equation 


(A-46)  ^  r 

^  J 


and  the  source  equation 


(A-47) 


c 


-  r} 


£ 


A' 


V 


Clearly  the  only  coupling  between  these  equations  Is  between 
Equations  A-43,  A-45,  and  A-46.  Substituting  Equation  A-46  into  A-43 
to  eliminate  yields 


Thus,  we  have  the  set  of  equations  given  In  Table  II  which  may  be 
solved  sequentially. 
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PABl  III:  Spironetcr  Equations  -  Table  III 

Both  pressure  and  temperature  In  the  spirometer  are  assumed 
constant.  The  temperature,  however.  Is  not  the  seme  as  that  In  the 
lungs  and  will  be  designated  T^.  From  the  equation  of  state  for  gas, 
we  have 

(A-49)  ^  C  ^  ) 

From  the  mass  balance  equation,  we  have  the  following  set  of 
equations : 


(A-50) 

c., 

(A-51) 

-- 

(A-52) 

'4.  -  k)  , 

and 

(A-53) 

Su^,  - 

■”  C^Urn  P^/z.  • 

These  equations  may  be  decoupled  by  combining  Equations  A-49  and  A-50 
to  yield 


(A-54)  K  ' 


r, 


- 


These  equations  may  now  be  solved  sequentially  as  llstsd  In  Table  III. 
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APPENDIX  B:  Stability  Test 

Within  the  context  of  finite  integration  a  differential  equa¬ 
tion  (  DE  )  la  convergent  if  the  solution  to  the  finite  difference 
equation  representing  it  approches  the  true  solution  to  the  DE  as  the 
independent  variable  step  size  approches  zero.  The  necessary  and 
sufficient  conditions  for  convergence  are  that  the  DE  be  both 
a)  conslstant  and  b)  stable.  Consistancy  indicates  that  the  finite 
difference  procedure  employed  does  infact  represent  the  DE  in  question 
and  not  some  other.  Consistancy  is  usually  assumed  and  is  only  exam¬ 
ined  if  the  equation  fails  to  converge  thou  it  proves  stable. 

Stability  proves  to  be  the  more  common  convergence  problem. 
Inherent  to  any  finite  difference  approximation  of  a  DE  is  some 
degree  of  error  including  truncation  error  and  round  off  error. 

There  exists  an  upper  limit  on  the  independent  variable  step  size 
below  which  these  errors  will  not  accumulate  over  a  number  of  inte¬ 
gration  steps,  in  which  case  the  equation  is  said  to  be  stable,  and 
above  which  they  will,  causing  instability.  It  should  be  noted  that 
stability  does  not  Indicate  that  the  error  will  be  small  -  only  whether 
or  not  it  will  be  amplified. 

The  equation  being  examined  is  the  ordinary  differential  equation 
(  ODE  )  for  pressure  in  the  airways,  which  in  its  most  general  form  is 

(B-1)  Pj  - 


^ [l  c  r  ^>7 
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No  convenient  technique  has  yet  been  developed  for  the  testing 
of  stability  of  the  ODE  though  several  methods  have  been  developed 
for  the  partial  DE  (  PDE  ).  It  now  becomes  a  problem  of  represent¬ 
ing  Equation  B-1  In  the  form  of  a  PDE. 

If  we  assume  for  simplicity  that  there  Is  but  one  Ch£unber  1 
and  one  Chamber  j  and  If  we  examine  Equation  B-1  at  time  t  using 


step  size  At  the  finite  difference  equation  would  be 

7  L  -V  - 


(B-2) 


^  ^ r  t  ^  c 

-p,  .  ^  f  -  /■',  c 


where  we  have  chosen  to  use  the  concentrations  at  the  current  time,  t  , 
In  the  approximation.  Recalling  that  concentration  is  a  function  of 
pressure  this  becomes 

NOW,  If  the  location  of  chamber  j  is  x  and  Chamber  1  Is 
(  X  -  AX  )  ,  Equation  B-3  Is  a  difference  analogue  for  the  PDE 

SP(tjx)  -  "•4'^ 

(B-4) 

A  stability  test  that  examines  the  PDE  was  developed  by  O'Brien, 
(46) 

et.  al. ,  and  has  since  been  treated  In  numerous  numerical  methods 
texts.  The  treatment  outline  advanced  by  Carnahan,  et.  al.,^^^^ 
was  used  to  examine  Equation  B-3.  One  assumes  that  a  separation  of 
variables  may  be  made  on  the  function  P(t,x)  and  that  the  space 
dependent  function  may  be  expanded  In  a  Fourier  expansion.  A  specific 
solution  may  then  be  chosen  from  the  Fourier  expansion  and  substituted 
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back  into  Che  PDE.  By  defining  Che  ampllflcaclon  facCor  f  repree- 
enced  as 


(B-5)  /  ^ 


/  {t 
T(i) 


where  T(c)  is  Che  Clme  dependenC  pare  of  P(C,x)i  and  where  Che 
condlclon  for  sCablllcy  Is 


(B-6)  f 

w 


<  / 


Che  upper  llmlC  of  Che  IndependenC  variable  ^C  nay  be  decemined. 

When  Chis  procedure  was  applied  Co  EquaClon  B-3,  Ic  was  found 
chac  i^\c  muse  be  less  Chan  zero,  which  was  Impossible,  indlcaclng 
ChaC  Chls  parclcular  flnlCe  difference  approxlmacion  was  absoluCely 
unsCable.  Ocher  approx Ima Cions  were  accempced  buC  none  could  be 
found  which  yielded  scablllcy  under  sacisfacCory  condlclons. 
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